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Abstract: We evaluate the contribution of chiral fermions in d = 2,4,6, chiral 
bosons, a chiral gravitino like theory in d — 2 and chiral gravitinos in d = 6 to all 
the leading parity odd transport coefficients at one loop. This is done by using finite 
temperature held theory to evaluate the relevant Kubo formulae. For chiral fermions 
and chiral bosons the relation between the parity odd transport coefficient and the 
microscopic anomalies including gravitational anomalies agree with that found by 
using the general methods of hydrodynamics and the argument involving the con¬ 
sistency of the Euclidean vacuum. For the gravitino like theory in d — 2 and chiral 
gravitinos in d — 6, we show that relation between the pure gravitational anomaly 
and parity odd transport breaks down. From the perturbative calculation we clearly 
identify the terms that contribute to the anomaly polynomial, but not to the trans¬ 
port coefficient for gravitinos. We also develop a simple method for evaluating the 
angular integrals in the one loop diagrams involved in the Kubo formulae. Finally we 
show that charge diffusion mode of an ideal 2 dimensional Weyl gas in the presence 
of a finite chemical potential acquires a speed, which is equal to half the speed of 
light. 
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1 Introduction 

The modifications of the macroscopic equations of hydrodynamics due to the presence 
of quantum anomalies of the underlying theory has been the focus of recent interest. 
It began with the observation of parity odd terms in the constitutive relation of 
the charge current in the holographic dual of J\f = 4 super Yang Mills at finite 
temperature and chemical potential [1-3]. This was then understood from general 
considerations using the equations of hydrodynamics, anomalous conservation laws 
in presence of background fields and the second law of thermodynamics [4] . These 
parity odd transport coefficients are non-dissipative . This fact was used to determine 
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the relation between microscopic anomalies and the macroscopic parity odd transport 
coefficient using a equilibrium partition function [5, 6]. 

Using general considerations from anomalies and the second law of thermo¬ 
dynamics or using the equilibrium partition function method it is not possible to 
determine the precise numerical constant that determine the relation between the 
parity odd transport coefficients and the gravitational anomalies or mixed gauge- 
gravitational anomalies of the microscopic theory. One reason for this is that both 
these methods rely on performing derivative expansions in velocities. Any relation 
between the gravitational anomalies or mixed gauge-gravitational anomalies would 
involve a jump in the derivative expansion. For example consider hydrodynamics in 
d — 2, the gravitational anomaly occurs at the 2nd derivative in the conservation 
law, however it affects the constitutive relations at the zeroth order in the derivative 
expansion [7]. There are three methods in literature which relates these anomalies 
to the parity odd transport coefficients. The first is the direct perturbative method 
of evaluating the transport coefficient of interest using the corresponding Kubo for¬ 
mula. This is usually done for free chiral fermions and then arguing that it is not 
renormalized, either using holography or in perturbation theory [8-11]. The second 
method is using the consistency of the Euclidean vacuum, developed by [7, 12]. The 
third method relies on performing one loop integration to obtain the theory on the 
spatial slice to relate the the anomalies to the Chern-Simons couplings and then 
arguing that these are one loop exact [11, 13]. The Chern-Simons couplings are in 
turn related to the transport coefficients 1 . To make this discussion concrete let us 
consider d = 2 and the the parity odd term in the constitutive relation for the stress 
tensor 2 given by 

= (e + P)u^u v - Prf v + A {2 \u^e yp u p + u v e w u p ). (1.1) 

Here u M is the velocity profile of the fluid. From general considerations of the equi¬ 
librium partition function and conservation laws [14, 15] 3 it can be shown that 
transport coefficient A® is of the form 

A^ = c 2d T 2 - c s /i 2 _, (1.2) 

where T is the temperature and //_ is the chiral chemical potential. The Euclidean 
partition method shows that c s is the coefficient of the charge current anomaly defined 
by the conservation equation 

d tl f = c s e^F, u , (1.3) 

F jiv is the background field strength. However the relation between c 2 d and the pure 
gravitational anomaly c g defined by 

d v T> a ' = FUf + c g ^V v R, (1.4) 

1 We thank Kristan Jensen and Zohar Komargodski for informing us about this approach. 

2 We have written down the stress tensor in the anomaly frame. 

3 See [16] for an alternative approach. 
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where R is the background curvature cannot be shown using methods which involve 
expansions in derivatives. This is because the coefficient c^d occurs at the zeroth 
order in derivative in the stress tensor. From (1.4) it is seen that pure gravitational 
anomaly is effective at the 2nd order in the derivative expansion. It was argued using 
an argument which involves the consistency of the Euclidean vacuum that 

hd = -8t r 2 c g . (1.5) 

A similar relation was obtained between the mixed gravitational anomaly and the 
parity odd transport coefficient in d — 4 determined by the two point function 
of the charge current and the stress tensor. In d — 4 the relation of the type in 
(1.5) has been verified by a direct perturbative calculation done in [8] as well as 
using the holographic dual of J\f = 4 Yang-Mills [9]. It was also shown within 
perturbation theory the relation in (1.5) is not renormalized for theories of chiral 
fermions which interact via Yukawa couplings. Similar equations relating parity 
odd transport and anomalies including gravitational anomalies were obtained for 
arbitrary even dimensions [7, 12, 17, 18] and summarized by the ‘replacement rule’. 
However the relation between gravitational anomalies and transport of the type given 
in equation (1.5) for dimensions other that of d = 4 has not yet been verified by a 
direct perturbative calculation. Further more it has been suspected that higher spin 
chiral fermions like gravitinos violate the relations of the type (1.5) [13, 19, 20] . 

Motivated by these questions we study the relationship between anomalies and 
parity odd transport coefficients using the direct approach of evaluating the respec¬ 
tive Kubo formulae. We consider theories of free chiral fermions in d — 2,4,6 and 
chiral bosons in d — 2 at finite temperature and chemical potential. We evaluate all 
the leading parity odd transport coefficients. This is done using the method devel¬ 
oped by [21] which involves evaluating the Feynman diagrams in finite temperature 
field theory and then performing the Matsubara sums. This approach is different 
form that is used in [22] to study transport properties of free chiral fermions in ar¬ 
bitrary even dimensions that uses Schwinger-Keldysh propagators We keep track of 
contributions that arise from gauge anomalies as well as pure gravitational anomalies 
and mixed gauge-gravitational anomalies. To address the issue of whether higher spin 
fermions obey the ‘replacement rule’ we study the contributions of a chiral gravitino 
like theory in d = 2 and chiral gravitinos in d = 6 to parity odd transport. 

We now briefly summarize the results of this work. The Feynman diagrams which 
contribute to the transport coefficients in d — 2 are considerably different from that 
in higher dimensions. This is because the spin connection is not present in the action 
of fermions in d = 2. From the perturbative calculation in d — 2 we show that chiral 
fermions as well as chiral bosons obey the relation in (1.5) predicted by the argument 
involving the consistency of the Euclidean vacuum. The transport coefficients of the 
chiral bosons is identical to that of chiral fermions which is expected from Bose- 
Ferrni duality in d = 2. For a higher spin fermion, that is a gravitino like theory, we 
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show that the the relation in (1.5) does not hold by an explicit perturbative analysis. 
In d — 4 we reproduce the results of [8, 21, 22] and show the perturbative results 
for the transport coefficients agree with the ‘replacement rule’. Since gravitinos are 
not charged, they do not contribute to the mixed anomaly in d — 4 and they do 
not contribute to any of the parity odd transport coefficients as well. In d = 6, 
for chiral fermions we reproduce the results of [22] for the contributions from the 
gauge anomaly to parity odd transport. We also keep track of contributions from 
gravitational anomalies. We show that contributions of chiral fermions to parity 
odd transport from the pure gravitational and the mixed anomaly agree with that 
predicted by the argument involving the consistency of the Euclidean vacuum. Next 
we consider gravitinos in d = 6 and show that its contribution to parity odd transport 
is not related to the pure gravitational anomaly and therefore does not agree with 
the prediction obtained from the replacement rule 

From out study of higher spin chiral fermions in d — 2 and d — 6 we can draw 
the following general lesson for their contribution to parity odd transport. These 
transport coefficients arise in correlators involving only stress energy tensor. We see 
that the reason that gravitinos do not obey the ‘replacement rule’ is because of the 
presence of the following term in their stress energy tensor. 

t 1 2 /2) “‘™ = _1 [ 8 ,(^ 7)A _^ r ) + &« o] • (i.6) 

This arises from the linearization of the Christoffel connection in the kinetic term 
for the gravitinos which is not present for fermions. Such a term gives rise to the 
interaction between the magnetic moment of the gravitino with the gravitational field 
[23]. It is this term that modifies the anomaly polynomial for gravitinos from that 
of a chiral spinor. However on evaluating Kubo formulae, at the zero frequency and 
zero momentum limit the contribution of the extra term in (1.6) vanishes. Thus the 
contribution of the gravitino to the transport coefficient is proportional to the chiral 
spinor. In fact it is just d— 1 times that of of the chiral spinor. The —1 accounts for 
the ghosts. We observe this phenomenon for the gravitino like theory in d — 2 and 
gravitinos in d — 6 while evaluating the Feynman diagrams involved in the Kubo 
formula, and note that this property is independent of dimensions. 

Apart form understanding the ‘replacement rule’ our explicit calculations also 
point to a technical simplification in evaluating the integrals involved in the one loop 
diagrams. We see that it is possible to take the zero frequency and zero momentum 
limit to obtain the transport coefficient before performing the angular integrations 
involved in the one loop diagrams. We develop a prescription for evaluating these 
integrals. This simplifies the calculation considerably compared to the ones in [8, 10, 
21]. In fact in [22] a guess is made for the angular integrals which was verified case 
by case. We will see that the angular integrals are elementary after taking the zero 
frequency and zero momentum limit. 


4 


Finally since parity odd transport coefficients appear at the zeroth order in the 
derivative expansion in d — 2 they can affect linearized hydrodynamic modes at 
the leading order. We show that the charge diffusion mode develops a velocity in 
presence of the parity odd transport coefficients. For the Weyl gas we see that this 
velocity is 1/2 times the speed of light. 

The organization of the paper is as follows. In section 2 we study parity odd 
transport in d — 2, we first obtain the relevant Kubo formula for the transport 
coefficient and then evaluate the contribution of chiral fermions, chiral bosons and a 
chiral gravitino like theory to the two leading transport coefficients. From the explicit 
calculation we will see how the ‘replacement rule’ breaks down for the case of gravitino 
like theory in d — 2. In section 3 we look at d = 4 and repeat the same analysis 
for chiral fermions. In section 4 we turn to d = 6 and study parity odd transport 
of both chiral fermions and chiral gravitinos. Again we see that chiral gravitinos 
do not obey the replacement rule. In section 5 we study linearized hydrodynamic 
modes which are modified by parity odd transport coefficients in d — 2. Appendix 
A contains details of moments of statistical distributions involved in obtaining the 
transport coefficients. Appendix B contains the details of the calculations of all the 
leading parity odd transport in d = 6. 

2 Anomalies and transport in d = 2 

In this section we first evaluate the 2 leading parity odd transport coefficients. They 
occur at the zeroth order in the derivative expansion in the velocities. To obtain the 
Kubo formulae for these transport coefficients, let us first write down the constitutive 
relations to the zeroth order in derivatives. The stress tensor and the current are 
given by 4 

= (e + P) - Prj^ + A {2 \u^e up u p + u u e^u p ), 
f = nuP + C (2) e^u v . (2.1) 

Here e, P, n refer to energy density, pressure and the charge density. A*- 2 -* and 
are the two parity odd transport coefficients which are of interest. The superscripts 
refer to the fact that these transport are for d = 2. Note that we have written the 
constitutive relations in the anomaly frame. We work in the metric with the signature 
(1, —1). To obtain the Kubo formula for these, consider the fluid perturbed from its 
rest frame tP = (l,v x ) with \v x \ « 1. We also perturb the metric and the gauge 
fields as follows 


9fu/ "F h pil , h tt , h tx 7^ 0) h X x 0, 

At A x — ct x . 

4 We work in the mostly negative signature through out this paper. 


( 2 . 2 ) 



Substituting this expansion in the constitutive relations in (2.1), we obtain the fol¬ 
lowing expressions to the linear order 


T tx = ev x + Pv x - Ph tx , 

f = -< (2 V + 


(2.3) 


Here we have equated the first order terms in the equations for the current and 
called it j* for convenience in notation. We have also used the fact that v x = 
—v x + h tx , to linear order. Solving for v x using the first equation and substituting in 
the constitutive relation for the current we obtain 


f = n- C <2 > 


+ PIhA 


+ i m h tI 


(2.4) 


This can be thought of as an Ward identity obeyed by the theory. We Fourier 
transform the equation and then differentiate with respect to h tx and finally set the 
perturbations to zero. Using the definition 


1 6 _ 

V9 $9hv 2 


we obtain 

{f(u,p)T tx (-uj,-p)) 


£(2) T tx (u,p)T tx (-uj , -p)) - P ^j _ ^ (2) 


( 2 . 6 ) 


Now we use the relation 


lim (T tx (c,p)T tx (-u,p)) = P, (2.7) 

which can be obtained by differentiating the first equation in (2.3) with respect to 
htx in (2.6). This leads us to the following Kubo formula 

- C (2) = lim (/(o;,p)T te (-u;,-p)) R . (2.8) 

p —40 

Similarly we consider the first order change in the stress tensor T tt from (2.1), which 
is given by 

T tt = 2\W(-v x + h tx ). (2.9) 

Again eliminating v x and differentiating the above Ward identity in Fourier space 
with respect to h tx we obtain 

= 1 lim (T“(u 1 ,p)T‘ I (-u,,-p)) R . (2.10) 

2 p -> o , uj->-o 

Before proceeding we make a few observations related to the Kubo formulae for 
the two leading order transport coefficients in (2.8) and (2.10). Both the two point 


6 






functions are obtained as the Fourier transform of the real time retarded correlator 
and the expectation value is taken over states held at finite temperature and chem¬ 
ical potential of the theory. These statements will hold true for all the correlators 
involving the various Kubo formulae studied in this paper. Also note that they do 
not involve any division by momenta, since these transport coefficients occur at the 
zeroth order in derivatives. We need to take the zero frequency limit first and then 
the zero momentum limit in (2.8) and (2.10). This is similar to the situation in higher 
dimensions where the static, zero frequency limit needs to be first taken and then the 
zero momentum limit [21], Note that these transport coefficients in d — 2 can also 
be obtained by considering one point functions of the stress tensor and the current. 
However we wish to illustrate the role of contact terms and also have a discussion 
similar to that seen in higher dimension and therefore we choose to examine two 
point functions. The transport coefficients of interest occur at the zeroth order in 
derivatives, they can in principle affect the hydrodynamic charge and sound modes 
5 at the zeroth order in the momentum expansions. We will study the consequences 
of these transport coefficients in section 5. 


2.1 Chiral fermions 

Our goal in this section is to evaluate transport coefficients (A 2 '* and A*' 2 -* given in 
(2.8) and (2.10) for a theory of free chiral fermions held at finite temperature and 
finite chemical potential. As we have discussed earlier, it is the Fourier transform 
of the real time retarded correlators which determine the transport coefficients. We 
obtain these correlators by first evaluating them in the Euclidean theory and then 
performing the necessary analytic continuation. Consider the Euclidean partition 
function S E of free chiral fermions, coupled to background metric and gauge held. 
The subscript E will always refer to Euclidean. It admits the following expansion 
around the unperturbed background given in (2.2) 


q _ c(0) I 

•JE - + 


A u + 


1 SS £ 

5 . 1 SS E 


2 y/gSg pa xfgSA, 


1 SS E 

\fg Sgpv 

)Apg pa + 


9p.v + 


1 5S h 


2 y/g6g pa Vgbg 

* 


)sy iv9po 




( 2 . 11 ) 


Note that due to minimal coupling, the gauge held A M occurs at most with a single 
power in the above expansion. The partition function of the theory is given by 


Z E = vp'vpexp (S E )- 


( 2 . 12 ) 


5 For relativistic fluids in 1 + 1 dimensions, these hydrodynamic modes were studied in [24], 
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Here i/j refers to the Fermion. Then the expectation value of the current and the 
stress tensor are defined by 


(f) E = 

(T^e = 


1 5 In Z F 


Vg SAp 

1 8\nZ f 


—0 


z 


( 0 ) 


\fg $9, 




-0 


z 


(0) 


V^V^e s E } 


V^Vxjje^ 


1 SSf 


J^5S e 

V9 ^9fj,u 


(2.13) 


9 lib 


Ruthin /—0 


where Z refers to the partition function in the absence of any perturbation. To un¬ 
clutter notations, from now all derivatives in the fields are understood to be evaluated 
at a^, h /iu = 0 . The relevant two point functions in the Euclidean theory are given 
by 


(fT^) E 


2 1 6 1 S\nZ E 

y/gSg pa 


(2.14) 


Substituting the expansion of the action given in (2.11) into the definition of the 
partition function we obtain 


( fT»°) E = -(f) E (T^) E + 


Z 


(o) 


T>A'D'ij)e“ E 


(o) 1 SS E 1 SS E 


z 


(o) 


1 6 


Vgfig pa Vg$A 


y/gSA^^/gdg pa 

1 SS E 


(2.15) 


Note the disconnected contribution of the first term cancels the disconnected contri¬ 
butions from the second term, in the first line of the above equation. The non-zero 
contribution of the correlator arises from the connected contribution of the second 
term and the from the last term. Since S E is a bilinear in the fermions, the third 
term is a contact term. We will see that such terms contribute crucially to the 
correlator and it is only after adding these terms, the ‘replacement rule’ for chiral 
fermions can be demonstrated. The structure of these terms in d — 2 is different 
from that in higher dimensions. On performing the same manipulations for correlator 
corresponding to A*- 2 ), we obtain 


(T^ T p a) E 


1 5 1 5 In Z E 

yfgbg pa y Vg $g (IV 


-(T^) E {T pa ) E + 


Z 


(o) 




e S E 


(o) 1 SS E 1 5Sz 


y/gtig (IV V9 5 9> 


pa 


Z 


( 0 ) 




e(0) 

e E 


1 8 , 1 SSp 


y/g&g pa Vgfig, 


)• 


(IV 


(2.16) 


Here again the disconnected terms cancel and we need to examine only the connected 
diagrams to evaluate the above correlator. 


























After evaluating the Euclidean two point functions, we perform the following 
analytic continuation to obtain the retarded correlator 

(j T (u n ,p)T TX {~u nr -p)) E = -(j\uj,p)T tx (-U, -p)) R \iw n ^po+ie- (2.17) 

Here c o n = 2nnT , where T is the temperature and n € Z refer to the Matsubara 
frequencies. Note that all external Matsubara frequencies are even multiples of nT 
and will be labelled by n. The overall negative sign results from converting the tensor 
indices containing the Euclidean time r to Minkowski time t using the replacement 
r —> it. Similarly we have 

(T xx (co n ,p)T TX (-co n , —p)) E = -i(T u (u;,p)T te (-u;,. -p)) R U^ Po+ie (2.18) 

Finally to obtain the transport coefficients, we use the Kubo formula in (2.8) and 

( 2 . 10 ). 

Let us now proceed to implement the above procedure in detail. The first step is 
to write down the expansion of the Euclidean action of chiral fermions in background 


gauge held and metric. The action is given by 

S E = [ drdx^gey^ a D,P^. (2.19) 

The gamma matrices convention, we choose to work with, are given by 

Y = ia 1 , Y = ia 2 , y c = — iYY — ~ (j3 - (2.20) 

The chiral projection operators are defined by 

P t =1(1=F7c). (2.21) 

Note that the Euclidean gamma matrices are anti-Hermitian, 7 C is Hermitian and P_ 
projects on to the upper component of spinor, while P + retains the lower component 
of the spinor. The gamma matrices obey the algebra 

{ 7 “,7 i } = 2 1 ,| t , 7 e=-V‘, (2.22) 

is the vierbien and the covariant derivative D M is defined by 

Dy. = df, + ^ajy Cd a cd + iAy, a cd = ^[YY} ( 2 -23) 


In two dimensions, the term involving the spin connection ui lic d vanishes and we do 
not consider it for the rest of this section. We will see in the subsequent sections that 
the contact terms in higher dimensions arise from linearizing the spin connection. In 
d — 2 the source of the contact terms arise from the linearization of the metric and 
the vierbein. 
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We proceed to expand the action (2.19) in terms of the following perturbations 
in the metric and the gauge held 


9tt 1 d h TT , 9rx 9 xt hrxi 9xx 1; (2.24) 

A t = A [ 0) + a T , A x = 0. 

Here At refers to the constant background chemical potential, which will be turned 
on subsequently. The inverse metric, to the second order in the perturbations is 
given by 


- \ — h - h 2 - h 2 -h - h h 

uv — I TT rr "'tx "'tx ' t rr' t ra \ D(h 6 

9 ~ 1 -h - h h - 1 - h 2 1 + ( 

it'TX Il'TT'l'TX - 1 - IL TX 


(2.25) 


Using the gauge in [23] the vierbein is given by 

f-l+^ + l(h 2 TX + h 2 TT ) _|_ hrrllTxX 




V 


h ’tx i 
~2~ ' 


TT ,b TX 


U 2 

1 I ll, TX , 
8 / 


+ 0(h 3 


Finally the inverse vierbein is given by 


(2.26) 


e f — 1 + + g (tf-x + ^rr) + 0(h 3 ), el — + -( h TX h TT ) + 0(h 3 ) 

e f = + g{h TX h TT ) + 0(h 3 ), el = 1 + -hl x + 0(h 3 ). (2.27) 


where the hatted variables refer to the hat space co-ordinates. We now substitute 
these expansions in (2.19) and obtain the following action, to the leading orders in 
the perturbations. 


S E 

s ( ‘> 


S‘e 


qW , C(2) 

1=1 E “T J 

J d 2 x^/g (x)Y^{x)^ x d x P^d>{x)\ 

— r Zl r(O ' > P_'0(x) — -^-^(x)^ x D t( -°' ) +7 T d x )P-d’{ x ) 


(2.28) 


3 h TT h T 


-{^(x)^ x D T ^ p_^[x) + ^ (x)Y d x P-d>(x)) 


ie / ddxyfg 




(a T '0 t (x)Y P-d>{x) + -d^a T '0' l '(x)7 :r P_'0(x) 


+ ... 

+ ■•• (2.29) 


where 

Z>(°) = d T + ieA?\ (2.30) 

is the covariant derivative in the background chemical potential. In (2.28) and (2.28), 
we have retained only the terms relevant to obtain the correlators of interest. Note 
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also, we have used the on-shell condition to simplify the coefficient of h TT . The 
components of the flat space stress tensor of interest is then given by 

T}[(x) = y YD T(0) P-it>, T}?{x) = + 7*(£ r,0) y-V (2.31) 

Similarly the time component of the flat space charge current is given by 

jji = ( 2 -32) 

Note that the quadratic terms involving the perturbations in (2.28) and (2.29) are 
responsible for contact terms in the two point functions of interest. These terms 
contribute in the correlators for d = 2 because the one needs only 2 gamma matrices 
to obtain a non-zero trace along with y c . As a simple cross check on these terms we 
considered the stress tensor in curved space and expanded it to first order in metric 
perturbations. The first derivative of the curved space stress tensor agrees with the 
second derivative of the action with respective to the metric obtained from 2.28 and 
2.29. 

The last ingredient we require to perform the computation of the two point 
function, is the propagator for the free chiral fermions. From the first term in the 
action given in (2.28), the Euclidean propagator at finite temperature and finite 
chemical potential in momentum space is given by 

(V>(^m,p)V’ t (kW,p / )) = S(u m ,p)6m,m'8(p ~ p')2n, (2.33) 

and 

S(uj m ,p) = -i[ 7 r (uj m - ei/j, - eifi 5 y c ) + YpV 1 , 

0 _ =i _ 

lUlm+e(p,-p. c )-p 

_ =i _ o 

tuj Tn +e(p+/j, c )+p 

Here c o m = (2 m + 1)tiT is the Fermionic Matsubara frequency at temperature T. 
The label m will refer to Matsubara frequencies which occur as odd multiples of ttT. 
We have chosen the background chemical potential to be 

A^ ] = iji + indci (2.35) 




where fi is the chemical potential corresponding to the vector 17(1) and /i c is the 
chemical potential corresponding to the axial 17(1) current 6 . Note that the trans¬ 
formation to the momentum space is given by 

x ) = ]}’£,[ ^(p) e - i{ “~ T * vx \ (2. 

—Y(p)e i(UmT+px) . 

2t r 

6 The signs for the chemical potentials are fixed using the thermodynamic definition dl g f f = j°- 
This allows us to identify the Euclidean gauge field to be Af = ifi. 
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This explains the relative sign between the frequency and the chemical potential in 
the propagator. 

£( 2 ) from chiral fermions 


P'+P 




Diagrams contributing to (7 2 ) 


Let us now evaluate the Euclidean correlator (j T T TX ). From (2.15) and (2.29), 
we see that there are two contributions. They are given by 


( f(p)T TX (-p )) 
hi = c r fl ( P )Tjn-p))c, 


A+B , (2.37) 

B = f / ^y 2 ^\p')l x P-4>{p')). 


Here, for convenience of notation, we have used p to refer to ( Po,Px ) = () where 
uj m is the corresponding Matsubara frequency. The integral in the second term B 
stands for 


[ d2p f— 

J (27t) 2 /3 J 2,71 


(2.38) 


Note that the term B arises from the second contribution in (2.15). The Fourier 
transform of the stress tensor and the current, in flat space, are given by 


3 }i(p ) = J -^y^iP -p)l T P~^{p'): (2-39) 

1 f d? z / 

Tjii-P) = 2 J + p "> ^ Tip ' x + 7 X ( i Po + + ^ c 7c))] P ~4>{p')- 


Substituting these expressions in (2.37) and performing the necessary Wick contrac¬ 
tions, we obtain for term A 

A = j J h TS (p' + P) ( 7 Vr + tWo + e (P + Pole))) S(p')P +}, ( 2 . 40 ) 

_ v- _e_ [ <¥x( _ P'x + ^ m + ep- \ 

^ 2/3 J 27r V (»w m + iwn + e/i_ - (p' x + p x )){iu m + e/r_ - pf x )) ’ 


where 


/i_ /i /i c . 


(2.41) 
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Note that there is an over all 27r/35(0) resulting from momentum conservation, which 
we have factored out. We perform the sum over Matsubara Fermionic frequencies by 
the usual trick of converting the sum to a contour integral. Here are the results for 
the sum 


y ~n [(iWm + e/i_ - p' x ))(iuj m + e/i_ + iuj n ~ p' x ~ Px )] 
m “ 

fip'x - e P-) - f(Px +Px- P-- iu n ) 


-1 


E 


1 


Px 

tUJm + G fl— 


j3 (iu m + e/i- - p' x ){iuJm + e/x- + iu n - p' x - p x ) 

Pxfjp’x - e/Q - (Px + P x - iu n )f(Px + p' x - iUn) _ 1 

2’ 


Px 

and / is the Fermi-Dirac distribution given by 

fix) = 1 


1 + e? x ‘ 


(2.42) 


(2.43) 


Note that the temperature, chemical potential independent constant 1/2 in the sec¬ 
ond sum can be ignored, since it just results in infinite constant. Substituting the 
results for these sums into the expression for A in (2.40), we obtain 


e f°° (¥xf 2p' x f(p' x -ep-) _ {p x + 2 p' x - iu n )f{p x + p' x - iu n ) \ ^ 

2 J —qo 27T y p x iio n p x iuj n J 


We can now change variables in the integration for the 2nd term and combine both 
the terms as ' 


ep x + iuj n 
2 Px ~ iUJ n 



)■ 


(2.45) 


Let us now examine the term H, again performing the required Wick contraction we 
obtain 


B = ~i \/ ^ T r(VS(p').P + ), (2.46) 

= jy [ dp * 1 

27T iU}m + e/i “ ~~ Gx ' 

After performing the Matsubara sum we obtain 

B = \ /_” - e?-). (2.47) 

'Though there is a divergence in the integrals when p' x —► — oo it can be shown that this infinite 
constant is independent of temperature or chemical potentials. The change of variables can be 
justified on careful treatment of the integrals when p' x < 0. 
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Adding the contributions in (2.45) and (2.47) the (j T T TX ) correlator is given by 


{f(p)T T ’(- p)> = ^ dp'J(p' x - e»J), (2.48) 

Px J — oo 

= e—[ ^(f(p> - e/i_) - f(p' x + e//_)), 

Px - iUn Jo 27T 

Px e 2 /i_ 

p x - iiOn 27r 

To obtain the second line of the above equation, we have again ignored an infinite 
constant, which is independent of temperature and chemical potential. The last line 
is obtained using the results of appendix A, equation (A.2), to perform the integral 
over the Fermi-Dirac distribution. We can obtain the retarded correlator using (2.17) 

(/(Po,Px)r te (-p 0 , ~Px))r = - —— e ~^- (2.49) 

Here ( Po,Px ) refer to the frequency and momentum in Minkowski space. We have 
ignored the ie that results from using (2.17), which is required when one needs to 
transform this correlator in Fourier space to a real time correlator. Finally using 
(2.8), the result for the transport coefficient is obtained by first taking the zero 
frequency limit p 0 —S > 0 and then the zero momentum limit p x —y 0. This results in 

C (2) = - lim (j t (po,Px)T te (-p 0 ,-Px))i?, (2.50) 

Px->0,po->-0 


2vr ' 

Examining (2.49), we see that it was crucial to take the zero frequency limit Erst 
and then the zero momentum limit. Note that both term A and the ‘contact’ term 
B in (2.37) contribute and they contribute equally in the zero frequency and zero 
momentum limit. Also note that in d — 2 term B contributed since one needed 
only 2 gamma matrices to saturate the trace. If once considers the contribution to 
this transport coefficient from the anti-chiral Fermions with chemical potential, then 
through a similar analysis one obtains 

A 2 ) — e U + (2 511 

^anti—chiral ^ 

where 


P-\- — P T Pc- (2.52) 

From (2.50) and (2.51), we see that the transport coefficient vanishes for a Dirac 
fermion in the absence of the chiral chemical potential p c . 
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Let us now compare the transport coefficient for chiral fermions in (2.50), with 
the coefficient of the U( 1) anomaly. The anomaly equation satisfied by the chiral 
current is given by 

d lt r = c s e^F, w , c s = -|^. (2-53) 

The value of c s for a Weyl Fermion in 2 dimensions can be read out at many places in 
the literature, see for instance equation (2.27) of [12]. We now compare the coefficient 
of the anomaly and the transport coefficient to see that they are related by 

C (2) = -2 c s n-. (2.54) 

which is the ‘replacement rule’ obtained in [7, 14, 15]. Observe that we have arrived 
at this relation between the transport coefficient and the microscopic anomaly from 
explicitly evaluating the diagrams that contribute to the Kubo formula. Note that 
the ‘contact’ terms gave rise to an equal contribution. 

For later purpose it is useful to parametrize the anomaly coefficients with the 
help of the anomaly polynomial. The anomalies of a theory in d dimensions, are 
encoded in a d + 2 form constructed out of gauge field strengths F and curvature 
tensor R ab . We define the curvature two-form R ab and the field strength two-form as 

Rab = ^Rabcddx c A dx d , F = ^ F ab dx a A dx b . (2.55) 

Furthermore, let us define a 2 /c-form which is a k-th degree polynomial in the curva¬ 
ture 2 -forms 

Rk= \R a a\ A £“3 A---^. (2.56) 

The anomaly polynomial for a single Weyl fermion in d = 2 takes the form [7] 

V d =i+i(F, R) = c s F A F + c g tr(R A R), (2.57) 


where, 


Co = —-—, c a =-. 

47t 967t 

y( 2 ) from chiral fermions 

We evaluate the correlator (r TT T TX ). Again from (2.16) and (2.28), we see that this 
correlator receives contributions from two terms which are given by 



(: T TT (p)T TX (-p )) = A + B, (2.59) 

A = (Tfi ip)Tfi{-p)) c , B = [7 x {Wt + (/i + ^c7c)) + YiPx} P-ip). 
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Figure 2: Diagrams contributing to 


To evaluate the term A we need 

T Jl ( P ) = J ^{p ~ P)Y{Wr + e(/i + /icT c))P-Hp)- ( 2 . 60 ) 

Inserting this and the expression for T TX (—p) from (2.39), and performing the Wick 
contractions using the propagator in (2.34), we obtain 

A(u n ,p x ) = 22 / ^r Tr (7 T [^(^n+ tu m ) + e(// + /vy c )] (2.61) 

x S(u> n + (Jm,Px + p' x )\iYp' x + 'ffam + e(/i + p c ^y c )]S (u m , p' x )P + ). 


Substituting the expression for the propagator, performing the trace and the Mat- 
subara sum results in 


j r°° dri 

lim A(0,p x ) = - -^p' x (f(p' x ~ ep-) + f(p' x + e/x_)), (2.62) 

p x ->-0 Z I n 27T 


t r 2 T 2 


" 8x 1 e2fl - + 


In the last line we have used the identity (A.4) to perform the integral. Evaluating 
B using the propagator given in (2.34), we obtain 


B = — 
8 f3 

f dPx , 

2i t 

n 

3 i 

r d Px 

" J, 

L ^ 

3 i 

{ 9 2 71 

87T 

1 e /i_ + 


2rp2 


(2.63) 


In the second line we have performed the resulting Matsubara sum, after which we 
have used (A.4) to perform the resulting integral. Finally adding the contributions 
from terms A and B results in 

hm(r-(0,p x )T-(0, -p x )) = ^ (e 2 p 2 _ + . (2.64) 

Using (2.18), the retarded correlator in Minkowski space is given by 

hm (T tt (0,p x )T tx (0, -p x )) = (e 2 p 2 _ + p-) ■ (2.65) 

Px~* o Z'K \ 3 J 
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Substituting this result in the expression (2.10), the transport coefficient of interest 
is given by 

A(2) ^ h ( eV - + ~ f) ■ < 2 - 66 > 

We now compare this transport coefficient with the coefficient of the gravitational 
anomaly of the Weyl Fermion, defined by the conservation law 


V„T^ = F£r + CgeTVvR, 



(2.67) 


The coefficient c g for a Weyl Fermion is obtained from (2.57). In [12], the transport 
coefficient A^ was parametrized as 


A (2) = c 2d T 2 - c s [i 2 _. 


( 2 . 68 ) 


Comparing with the result from evaluation of the Kubo formula given in (2.66), we 
see 


c, = 


(2.69) 


e„ 7r 
’47t’ Cm_ 12' 

The value for c s agrees with the chiral anomaly of a Weyl fermion in 2 dimensions 
as in (2.53). The value for c 2 d obeys the relation 


C 2 d = -87 T 2 Cg. 


(2.70) 


as argued in [12] using the consistency of the Euclidean vacuum. We have arrived at 
these relations by explicitly evaluating the diagrams which contribute to the Kubo 
Formula. Note that for both transport coefficients, there was a contribution from 
the ‘contact’ term which was important in obtaining these relations. 


2.2 Chiral bosons 

In 2 dimensions a massless chiral Weyl fermion is dual to a chiral boson. Therefore 
we expect the contributions of a chiral boson to the transport coefficient and 
A^ 2 ), to be identical to that of the chiral fermions studied in section 2.1. In this 
section we will show that it is indeed the case. Though the result is expected on 
general grounds, performing the calculation in detail will serve as a cross check of 
the analysis in section 2.1. Also as far as we are aware transport coefficients have 
never been evaluated directly for chiral bosons. With these motivations we proceed 
with the analysis of transport in this system. 

Chiral bosons in the Euclidean theory obey the constraint 8 

I1U 

d^(j) = (2.71) 


e TX = e TX = 1 and is anti symmetric in its indices. 
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Note that this is the constraint which in flat space reduces to d s (j) = 0 on using 
the change of variables to holomorphic and anti-holomorphic coordinates given in 
(2.109). The constraint ensures that the bosons depend only z, this property is is 
identical to the Weyl fermions considered in section 2.1. To write down the stress 
tensor and the current of the chiral boson theory, we impose the constraint of (2.71) 
in the stress tensor of the boson and the current following [23]. They are given by 

, (2.72) 

If \ 

f = (tr + i—a^j, 

where e& is the charge of the boson. The normalization of stress tensor is consistent 
with that used in conformal field theory in two dimensions. We have kept the metric 
in these expressions for the stress tensor arbitrary, to enable us to evaluate the 
relevant Kubo formula. The partition function of the theory is identical to that of 
chiral Weyl fermions. The frequencies for the bosons in the Euclidean theory are 
even multiples of ttT and are given by 

oj n = 2nnT 1 n G Z. (2.73) 


rjiflU _ _ 

~8~ 




d^ + i—dp, 

V9 


0 <9 


+ i 


The thermal two point function for the bosons is given by 

. 2 

S B {u n ,p) = (0(cu n ,p)0(-o; n /, -p )) = -—- -2n(35 n y5(p - p ). (2.74) 

[iu n y -p z 

An important fact to be remembered while performing Wick contractions is that 
the one point function of the charge density in the thermal state is non zero in the 
presence of the chiral chemical potential and is given by 

(j%x)) = e^. (2.75) 

This can be derived from the partition function of the fermion theory which coincides 
with that of chiral bosons 9 . One can show that the expectation value of the charge 
density is given by (2.75) also directly in the bosonic theory by adding the term /i_j° 
to the Hamiltonian and evaluating the partition function. Converting the expectation 
value of the charge to the Euclidean theory we obtain 

U t (t,x)) = ie 2 ^. (2.76) 

Let us convert this expectation value to a rule for determining one point functions of 
the operator (j) in momentum space. Substituting the expression for j T from (2.72) 
and performing the Fourier transform to momentum space we obtain 

2 

{(iu n + p x ){<l>(u n ,p)) = 27 T/35 n , 0 5(p) — — . (2.77) 

e b 7t 

9 We review the thermodynamics of the Weyl gas section 5. 
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We will now determine the relation of the charge to the charge of the chiral fermion 
e by demanding that the thermal expectation value of the stress-tensor T TT agrees 
with that of the Weyl fermions. The stress tensor in momentum space is given by 


Tjl ( UJ n , p) = J l iUJ n' ~ iUJ n) + ( p' - p)] l iuJ n' + P') ( 2 ' 78 ) 

X0(cU n /,p / )(/)( U n -u n ',p-p'), 

The thermal expectation value can be written as, 


(T rr ) = A + B. 


(2.79) 


= s ~ iuJn ) + ~ p ^ [' iUn1 + ( 2 - 8 °) 

n' 

x(0(uv,p')0( uj n -u n ',p-p')), 

= J ^ iuJn ' ~ iuJn "> + $ ~ p )] \- iuJn ’ + p/ ] 

n' 

x(0(ov,p'))(0( u n ~u n ',p-p')). 

We now evaluate term A by substituting the propagator from (2.74), this yields 

a _ 1[ dp ’ [ iu n '+p'} 2 

4 P J 2vr (iav) 2 - (p') 2 

—hljMW 


A 


B 


Here we have factored out the 2n5 nt0 5{p) which results from momentum conservation. 
b(p x ) denotes the Bose-Einstein distribution. To obtain the second equality we have 
used the following Matsubara sum for even integer frequencies 



1 iC0n' 

P ( iuj w - p' x ) 


Kp'x) 


-p'xKp'x ) 
1 

e Pp'x - 1 ’ 


(2.81) 


On neglecting an infinite constant independent of temperature we obtain 


A = ~— dp p' b(p). 


(2.82) 


Using the Erst moment of the Bose-Einstein distribution given in ( (A.7)) we obtain 

7 tT 2 


A = 


12 


(2.83) 
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Let us now examine the term B in (2.80) Substituting (2.77) for the expectation 
values of the currents we obtain 


e g 7T 


(2.84) 


Here again we have factored out the delta function due to momentum conservation. 
Then putting terms A and B together we obtain 

(2.85) 


1 tt 2 T 2 e A u? 

( T TT\ = 

[ ' 47t 3 e?27r' 


Comparing this to the energy density of chiral fermions given in (5.13) 10 we see that 
for Bose-Fcrmi duality we must define the charge of the chiral boson to be given by 

e 


e b = 




( 2 . 86 ) 


Now that we have defined the theory we can proceed to evaluate the transport 
coefficients. 


£(2) f r om chiral bosons 



Figure 3: Diagrams contributing to Ccwraiboeon 


To obtain the possible contact terms, we expand the expression for for the j T 
component of the current, to the linear order in the metric. The relevant terms are 
given by 

j T = yMr + id x )<j> - ^pL(-<9 T + id x )4> + 0(h, h 2 ). (2.87) 

Then the is obtained by examining the following correlator 

+ < 2 - 88 > 

10 Note that the overall negative sign in ( T TT ) is due to the fact we are in Euclidean space. 
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Let us now evaluate the correlator A by performing appropriate Wick’s contractions. 
In momentum space we have 


jfi(u n ,p) 


T^{-u n ,-p) 


— (iu n + p)(p(uj n ,p), (2.89) 

i J + P ) ^ n ' + iuJn ^ + (P + 

n' 

X0(ca n /,p , )0(-( w n + w n /), -(p + p'))- 


For the present correlator, due to the structure of the Wick contraction and the 
momentum and conservation, there is no sum over internal momenta. Evaluating 
the two possible Wick contractions we obtain 


e 2 p_ iu n + p 
Air p — ico n 


(2.90) 


Here we have used the expectation value (2.77). The diagrams contributing to A are 
shown on the LHS of figure (3). The contact term B is given by 


B = — (iu n +p)(0(cn n ,p)) 


e 2 /r_ 
47T ’ 


(2.91) 


where again we have used (2.77). The diagram contributing to B is shown on the 
RHS of figure (3). Note that as before we have suppressed the 27 t/ 3<5(0) terms in A 
and B. Adding the two contributions we obtain 


C f(p)T TX (-p )) 


p e 2 /i_ 
p — iu n 27r 


(2.92) 


This is identical to the expression obtained in the fermion language given in (2.48). 
Therefore taking the zero frequency limit first and then the zero momentum limit 
we see the contribution to the transport coefficient from the chiral bosons is 
identical to that for the fermions and is given by 


A2) = eV 

chiralboson r* 


(2.93) 


Note that the contact terms contributed equally just an in the case of the chiral Weyl 
fermion. 


y( 2 ) from chiral bosons 

Let us expand the T TT component of the stress tensor to the linear order in the 
metric. The relevant terms to obtain the contribution of the contact terms are given 
by 

T TT = —-(1 + ih TX )(-d T + id x )(t)(-d T + id x )<f> + 0(h , h 2 ). (2.94) 

O 
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p'+p 



The correlator from which A*- 2 -* is extracted is given by 


1 5 f 6\nZ 

Jgbh 

TX ip) \5h TT (~P) 


A = (TJ[{p)TJi ( —p)), B = 


= A + B , 

, ST”(p) 

'Sh TX (p ) ' 


(2.95) 


In term A and B we need to evaluate all the contributing connected diagrams. The 
components of the flat space stress tensor in momentum space are given by 

(w„, p) = l 2^ ~ iUn ) + ^ _ p )] + ( 2 - 96 ) 

“ m J 

X(j)(uj m ,p')(j)( UJ n - U m ,p- p'), 

x - ^ /* fyv/ 

Tfl{~U n , ~P) = / 2^( iUm + P ^ ^ iUm + iw **) + (P + P')] 

X(j)(u m ,p')(j)(-(uJ n + cu m ),-ip + p')). 


Let us proceed to evaluate the term A. This quantity has two contributions due to 
the fact that the current has an expectation value in the thermal state. 


A 

A! 

A 2 


A\ + A 2 , 
2 i 


(2.97) 


I 7^(^n' +p') 2 {iu' n - UJn + P ~ p) 2 S B (u n ', p) S B {u n ~ U n >,P~p'), 

n I 


j' . X! / [{iWn>+p') 2 {i(J n -Un+lf -p) 2 S B (Un',rf) 

" n' 

x{(j){u n ^u' w p-p'))] (0(w n / -UniP' -p). 


These two contributions are represented by the diagrams on the left hand side of 
figure 4. Note that there are 2 possible Wick contractions which contribute equally 
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giving rise to the factor 2 in A\ and there are 4 possible Wick contractions which 
contribute equally giving rise to the factor 4 in A 2 . Substituting the expression for 
the bosonic propagator Sb from (2.74) in to A\ we obtain 


A\ — 



dp' (icy + p') (icy - iuj n +p' -p) 

2 vr (icy - p')((iu> n ' - iu n ) -(p'-p))' 


(2.98) 


We can now perform the Matsubara sum and take the zero frequency limit first and 
then zero momentum limit to obtain 


lim A 1 (p 0 ,p) = l - f yV&(?/), 
p—2 J 2ir 

iirT 2 

~ 12 


(2.99) 


From the structure of the Wick contractions in A 2 we see that it does not involve 
any Matsubara sum or integral. We need to use (2.77) to evaluate the expectation 
value of the currents and the relation between the charges in (2.86). This leads to 
the following value for A 2 



( 2 . 100 ) 


The contribution from the contact term B is given by 


ST”(p) 

Sh TX (p ) 


i 

47T 


(e 2 /i- + 


3 ’ 


( 2 . 101 ) 


This is easy to see from (2.94) since the contact term is equal to —i(T TT ). The 
diagrams contributing to B are shown on the RHS of figure (4). Considering the 
contributions both from the term A and B we obtain 

lim^yto.pjpr(0 - p)> = ^(eP_ + 7T). (2.102) 


This is identical to that the contribution of the chiral fermions for the corresponding 
correlator as seen in (2.64). Therefore by the same logic we obtain the contribution 
of chiral bosons to A 1 - 2 ) to be given by 

■^chiralbosons = ^ -)' (2.103) 

It is indeed interesting to see how the different Feynman diagrams in the chiral boson 
theory which obeys the Bose-Einstein distribution organize themselves to give the 
same transport coefficients as that of the chiral fermion theory. 
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2.3 Chiral gravitino like system 


It has been suspected that fermions with spin greater than 1/2, for instance graviti- 
nos, do not obey the replacement rule [19, 20]. That is the relationship between parity 
odd transport and the microscopic anomaly breaks down. It will be instructive to see 
how this explicitly occurs in two dimensions, by performing a perturbative evaluation 
of the respective Kubo formula. There are no physical propagating gravitinos in two 
dimensions. However the following system mimics the features we require. Consider 
the Lagrangian 


c _ ca i I cs 2 
de — o e ' J E i 

S S e = J drdx^eylrD.P-V, (2-104) 

where 

D^ p = + -uj^cd(y cd V + T (2.105) 

This Lagrangian in (2.104) is the gauged fixed Lagrangian of the gravitino in higher 
dimensions. The action S 9 ^ is the action of the corresponding ghosts. The con¬ 
tribution of the all the ghosts to the gravitational anomaly can be accounted by 
subtracting the contribution of a chiral spinor [23]. In two dimensions the above 
Lagrangian would not have any physical propagating degrees of freedom. We will 
see that in fact the action S f- 1 has central charge c = —22. However this gravitino 
like system will provide a simple illustrative example of how the replacement rule 
can break down 11 . 

To proceed let us evaluate the central charge of this system. There are additional 
terms in the stress tensor due to the presence of the Christoffel symbol in the covariant 
derivative given in (2.105). The linearization of the Christoffel symbol results in 

r p« = ^rf’idphtra + dahap - d„h ap ) + 0(h 3 ). (2.106) 

The stress tensor of the system is given by 

Tpv = T 9 l + T 9 l, (2.107) 

T pl = + lud^P-Vix) 

T fv = pP-M x ) ~ <Wi(x)7 vP~V{x) + (|iHr)). (2.108) 

Note that the total derivative term in the stress tensor arises from the contribution 
due to the linearization of the Christoffel symbol. While evaluating the stress tensor, 

11 We thank R. Loganayagam and E. Witten for discussions which clarified our understanding of 
this point. 
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we have used the equations of motion. To make it convenient to determine the central 
charge we go over to holomorphic coordinates defined by 


z = — ir + x, z = ir + x. (2.109) 

Then the only non-trivial component of the stress tensor is given by 

Tzz(z) = ^ (^*3^ - < 9 ^*^) - ( 2 . 110 ) 

Here b> /t is the chiral component of the gravitino, e is the anti-symmetric tensor 
defined with e TX = 1. We can now easily evaluate the central charge with the following 
0PE 

(2.H1) 

p z — w 

On performing this, it is easy to see that there are additional contributions for the 
term proportional to the total derivative in (2.110). The central charge is given by 

4 1 = 2 - 24 = -22. (2.112) 

The 2 results from the first term in (2.110), while the —24 results from the total 
derivative. Now to obtain the complete contribution, we need to take into the con¬ 
tribution from the ghosts, which is equivalent to subtracting the contribution from 
a complex chiral fermion. Therefore we obtain 

c R = 4 1 + 4 2 = -22 - 1 = -23. (2.113) 


The reason that this central charge is negative is due to the fact that the gravitino 
like system in 2 dimensions is not physical. The contribution of this gravitino like 
system to the gravitational anomaly is given by 


, gravitino _ 


23 


~ 9 96t t Cr 96tt' (2.114) 

Here we have used the relationship between the central charge and the coefficient of 
gravitational anomaly form [7]. This contribution to the gravitational anomaly of 
the gravitino like system coincides with what has been considered as the contribution 
of the ‘gravitino’ in [19]. 


y( 2 ) from gravitino like system 

We now evaluate the contribution to the transport coefficient A® from the grav¬ 
itino like system. The components of the classical stress tensor for this system in 
momentum space can be read out from (2.107), they are given by 

Tflip) = T T ,; m + TJf m , (2.115) 

T n x) Cp) = / -pmipr)p-rw), 

77 ,( 2 ) (?) = ip, J 477 [«■'"(?' - ph T p~r(p')) - (P r (p' - P )YP~r(p'))\ . 
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( 2 . 116 ) 


rpTX ( \ _ rprx j rpTX 

1 fi (-P) ~ 1 fl( 1) + J /z(2)> 

T m(-P ) = | (?' + p)(7V, + 7>r)P-f (?'))). 

T fli2)(-p) = J Jr^y [(^ a (p' +P)YP-Vip') +^ a (p' +p)j X P-^ T (p')) 

- (^ T (p + p)^ x (p) + •0 ta: (p / +p)7 T -P-'0 CT (p / ))] • 

Here the terms in the stress tensor, which arises from the total derivative, have 
been separated and their contribution is labeled by the subscript (2). The two point 
function of the classical part stress tensor is given by 

(Tfl(p)T}f(-p)) = (Tj[ m (p)T]? m (-p)) + (T/7(p)ry 2) (-p)). (2.117) 

It can be easily seen from the structure of the Wick contractions, all the remaining 
cross terms vanish. Now the contribution from (Tj^ j^ (p)Tjfi j T—p)) is twice that of 
a single chiral fermion. We will show that the contribution from (TJ^(p)TJ^ 2 )(—p)) 
vanishes. We evaluate this correlator by performing Wick contractions using the 
propagator 

(i> t i(^m,p)'fpl{^m',p') = - p')rJnvSm,m' ( ^ ' Um ~ P j • (2.118) 

\ iUm+P J 

Note that we have considered the gravitino like theory at finite temperature but 
not at finite chemical potential. This is sufficient to pick out the pure gravitational 
contribution to the transport coefficient A* 2 ). After performing the Matsubara sums 
involved we obtain 

( T f[{2){u m ,p)Tj? {2) (-uj, - p)) = W ^' n+ _ P n [ Try (fip) ~ f(p'+p)), (2.119) 

ICUm P J —oo ^ 

. P 2 iu m +p 

2tt ico m — p 

Therefore we see that in the zero momentum limit the contribution from the Christof- 
fel symbol vanishes in the zero momentum limit. 

n lim A T m( u m,p)Tj? {2) (-u t -p)) = 0. (2.120) 

For later analysis it is important to note that even if we did not evaluate the integral 
in (2.119), as long as one is interested in the finite piece of the integral, taking the zero 
frequency limit and then the zero momentum limit ensures that the this contribution 
vanishes. Thus we conclude that the contribution to the transport coefficient A® 
of the total derivative term in the classical stress tensor (2.107), which arises due 
to the presence of the Christoffel symbol, vanishes. Let us now examine possible 
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contact terms which arise from expanding the Christoffel symbol. On performing 
the expansion of the vierbein, metric as done in (2.25), (2.26) and (2.27) as well as 
the term involving the Christoffel symbol in the action (2.104), the relevant contact 
terms arise from 


o(i) _ 

°Christoffell — 


clrdx ( h TX d x h TT (^l r fij) x ) + 


h a h 

,l TX^T , l TT / , f XI \ 

-x- Wrl Wt) 


h r) h 

,b TX w X ,lj TT/ , + T I \ 

---(Vv7 Vv) 


( 2 . 121 ) 


Note that the origin of all these terms is from the expansion of the Christoffel symbol 
and therefore they involve a derivative on the metric. This is unlike the other terms 
in the second line of given in (2.28), which involve a derivative on the fermion. 
The presence of the derivative on the metric leads to a momentum factor outside 
the Fermi-Dirac integral after performing the Wick contraction and the Matsubara 
sum. Therefore again, on taking the zero frequency and zero momentum limit, the 
contributions of these terms vanish. 

The above analysis shows that the terms arising from the Christoffel symbol in 
the action for the gravitino like system in (2.104), do not contribute to the transport 
coefficient A^ 2 b Thus the contribution of this system to this transport coefficient is 
identical to two chiral Weyl fermions. We then have to subtract the ghost contribu¬ 
tion which is equal to a single chiral Weyl fermion. The net result is that A* 2 ^ for a 
gravitino like system is given by 


A (2 > = 


7 tT 2 

~Y2 


( 2 . 122 ) 


Therefore we have 


~ gravitino 
C 2d 


7r 

12 ' 


(2.123) 


Comparing the value of c ^s ravltmo for the system in (2.114), we see that this sys¬ 
tem does not satisfy the replacement rule obtained by consistency of the Euclidean 
vacuum since c 2 ® ravltmo ^ —87r 2 c 9 gravltino . 

At this point it is important to make the following observations. The terms 
proportional to the Christoffel symbol contribute to the anomaly coefficient. However 
for the transport coefficient the contributions of this vanished on taking the zero 
momentum limit. It is for this reason that for the purpose of evaluating the transport 
coefficient the system along with the ghosts behaved like a single chiral Weyl fermion. 
Though the system we analyzed in two dimensions possibly does not arise naturally 
in any physical system since the central charge of this system is negative, it illustrates 
the mechanism of how for higher spin chiral fermions the transport coefficient is not 
constrained by the microscopic anomaly. We will return to this phenomenon again in 
6 d dimensions, where chiral gravitinos are physical and they contribute to the pure 
gravitational anomaly. 
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3 Anomalies and transport in d = 4 


Let us begin by parametrizing the constitutive for the stress tensor and the charge 
current. We will restrict our attention to terms occurring in first order in the deriva¬ 
tive expansion in the anomaly frame. In d — 4 the constitutive relations are given 
by 

T pv = (e + P) u»u v - PrT + AS A \u^e vpaP u p F aP + u^e^u p F aP ), 
+\¥\u»e upa0 u p d a up + u v e ppap u p d a up), 
f = nu p + Ci (4) \e ppap u p F a p + d 4) e ppa \d a u p . (3.1) 

where, e,p, n, /i, Ai, A 2 , Ci, C 2 refer to the energy density, pressure, charge density, 
chemical potential and the four anomalous transport coefficients respectively. The 
metric, we work in, has the mostly negative signature g pu = (1, —1, —1, —1). Let us 
derive the Kubo formulae for the parity odd transport coefficients, by considering 
perturbations of the fluid from the rest frame. Let the velocity of fluid be given by 
u p = (1,0, 0) with \v y \ « 1, we also have such that u p u p = 1 + 0(v 2 ). The 

gauge held, velocity and metric perturbations are given by, 


A p = (0, 0, a y , 0,) u p = (0, 0, -v y + h ty , 0), 

hfiui hty hty 7 ^ 0. (3.2) 


All perturbations are assumed to depend on only the z direction and we can work 
with the Fourier modes in this direction. Expanding the constitutive relation to 
the linear order and interpreting it as an Ward identity, just as it was done in the 
case of d — 2, we obtain the following Kubo formula for the anomalous transport 
coefficients. 


A< 41 = 
A< 4, = 

c! 4, = 

d 4, = 


lim 

p z -s>0,po->0 


lim 

Pz- 5 > 0 ,po-aO 


lim 

Pz ^0,po ^0 


lim 

Pz AO,po ^0 


— (T tx (p 0 ,Pz)j y (-po, ~Pz))r, 

^Pz 

— (T tx (po,P z )T ty (-p 0y -p z ))R 
IPz 

— {f{P0,Pz)j y (-P0, ~Pz))r, 

Wz 

-(j X (P0,Pz)T ty (-p 0 , Pz)) R- 
^Pz 


(3.3) 


The expectation values in all these correlators refer to real time retarded correlators 
at finite temperature and chemical potential. As discussed earlier, we have to take the 
zero frequency limit first and then the zero momentum limit. All of these correlators 
involve division by a power of momentum, since these transport coefficients occur at 
the first order in the derivative expansion. Finally note that the transport coefficients 



A ( | 4) and rf 4) are identical and therefore it is sufficient to evaluate only 3 of the Kubo 
formulae in (3.3). As done for the case of d = 2 to obtain the Minkowski real time 
retarded correlators, we will evaluate the corresponding correlator in the Euclidean 
theory and then analytically continue to obtain the transport coefficients. 

Before proceeding, let us recall previous evaluations of these transport coefficients 
in the literature. For a theory of Weyl fermions, the coefficient was first evaluated 
in [21], the coefficient Cl 1 * as well as was evaluated in [10]. As we will see, the 
coefficient A, 4 ^ is sensitive to the mixed gravitational anomaly in the theory and 
this was obtained in [8]. The evaluation of the coefficient Ag was begun in [10], 
and all the diagrams including the contact terms which contribute to this coefficient 
was finally done in [22], The fact that the mixed anomaly contributes to A, 4 \ was 
mentioned in [10]. However the complete evaluation of the contribution of the mixed 
anomaly coefficient including the contact terms has not been done. The work in [22] 
was interested only in the contribution of the global gauge anomaly. We will complete 
this small gap in literature, we will also develop a simple method of evaluating the 
resulting angular integrals. Our calculation will also serve as a cross check of the 
methods in [22] which are different from that in this paper. We also show that 
we can take the zero frequency and zero momentum limit Erst and then perform 
the resulting angular integrals. This considerably simplifies the evaluation of these 
integrals. In all the previous work, the angular integrals were first performed, in 
fact [22] provides a formula which was guessed for angular integrals which occur in 
arbitrary dimensions. 


3.1 Chiral Fermions 


Let us proceed to evaluate the corresponding correlators first, in the Euclidean theory 
of free chiral Weyl fermions in d = 4. Following the set up in section (2), the 
expression for the Euclidean (j p j u )e is given by 


{ff)E,cl = -(f) E {j U )E + 


z 


(o) 


v^v^ 5SeSSe 


da M 5a h 


(3.4) 


Z 


1 , e(0) 5 2 Se 

T>^T>t[je SE 


(o) 


5a^8a h 


while the two point functions (j p j u )E and (j p T up ) e are given in (2.15) and (2.16) 
respectively. For convenience, let us define the following coefficients, which are ob¬ 
tained from taking limits of the Euclidean correlators 

S 4) = lim n — {j x (Q,Pz)j v (0,-Pz))E, (3.5) 

p z ->o, xp z 

c < 41 = lim r—(j x (0,p z )T Ty (0, — Pz))e, 

p z ->0 ip z 

A< 41 = lim —{T”(0,pJT'»(0,- P «)) E . 

Pz-> 0 ip z 
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To obtain the transport coefficients from the above definitions, we use the relation be¬ 
tween Minkowski real time retarded correlators and the Euclidean correlators which 
are given by 

(j X (uj n ,p)j y (-UJn,-p))E = ( j X (po,p)j y (-po,-p))R\iu^ Po+ ie , 
(j x (u n ,p)T Ty (-uj n , -p)) E = i(j x (p 0 ,p)T ty (-p 0 , -p)} R \i^p 0+ ie, 

(T TX (u n ,p)T Ty (-uj n , -p)) E = -(T tx (po,p)T ty (-p 0 ,-p)) R \^ po+ie . 

(3.6) 

where, u> n = 2nnT, are the integer quantized Matsubara frequencies. Using (3.6) 
and the definitions of the transport coefficients in (3.3) and (3.5) we obtain 

ci 4) = d 4) . ci 4) = id 4) . a' 4 » = -a< 4 ». (3.7) 

For Weyl fermions in d — 4, the partition function Z E is given by 

Z E = j Vi/}lpip E exp{S E ), (3.8) 

where, 

S E = j drd 3 x^ge^'y a DpP_'ilj. (3.9) 

The covariant derivative is defined as 

Dpif} = d p ip + ^u pab (T ab '^ + ieApip, (3.10) 

where, u) pa b is the spin connection and a ab = ^[7 a , 7&]. Unlike d = 2, the contributions 
from the term proportional to the spin connection does not vanish for d > 2. We 
work in Euclidean space where the gamma matrices satisfy, 

{ 7 “, T ‘} = 27 ? g>, — —6 ab , (3.11) 

and 

e±=i(l±75>. (3.12) 

Note raising and lowering of flat space indices are done by . We now perturb the 
action in (3.9) from flat space and constant chemical potential, by considering the 
expansion with the following choice of metric and gauge field perturbations, 

9pv — —dpi; + hp U , Ap = + a,p, (3.13) 

where At refers to the background chemical potential. For this metric perturba¬ 
tion following [23], we can work with the following vierbein to the linear order in 
fluctuations 

Gap — —dap H—+ 0(h 2 ). (3-14) 
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Essentially we are working in a gauge which ensures that the vierbein is symmetric. 
Expanding the action to the quadratic order in the fluctuations we obtain 


S E = S$ ) +S$ ) + S%\ (3.15) 

4 0) = f d^iYDr+fdJP-'l/jE, 


S ( E ] 

q( 2) 

o E 


— J + Y^)P-^e + ie 

J 16 


d^x^YA^P-^Ei 


(3.16) 


Here coefficient of the term linear in the metric is the stress tensor in flat space and 
the coefficient of the term linear in the gauge field is the charge current. The second 
order terms in the the metric, given in the last term of (3.15), arises from expanding 
of the spin connection, the repeated indices in the subscript refer to summation. This 
term contributes to the mixed anomaly coefficient [23]. We will see that that this 
term gives rise to a contact term and it is important in the evaluation of the coefficient 
y( 4 ) va s t anc i s f or totally anti-symmetric combination of gamma matrices which 
is given by 

Y ua = i (7^7“ - 'fY'f + (cyclic) ). (3.17) 

Using the action (3.15), we can obtain the stress tensor and charge current which 
are given by 

7TM = \^ E (Yd- + YWP-I’e, (3.18) 

jfi = -ie^ E YP-^E- (3.19) 

We now take the background chemical potential to be 

= i/i + i/vys- (3.20) 

The propagator derived from the action (3.15) is given by 

—i 

Y [w m - ien - ie/vy 5 ) + 7 % ’ 

A t(i^s,q)P s r y^t, (3.21) 

s,t=± 

?± = (-*>±?), ^ 9 = lei (3-22) 

cog = co m — ie/j, — siefi c , ji± = /1 ± /i c . 


Se (q) = 


where 


P + = 


A ±(qo,q) = 


(! ± 75 ) 


2 

1 


qo T Eq 
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p'+p 



Figure 5: Diagram contributing to (j 4) 


Evaluation of 

From the Kubo formula in (3.5), the transport coefficient is given by 

= lim ——{j x (0,Pz)j y (0 1 —p z ))- 

Pz^f o ip z 

(3.23) 

For all calculations of correlators in the rest of the paper, we will set the frequency to 
zero at the beginning. The correlator is defined using (3.4) and from the expansion of 
the action given in (3.15), we see that there are no contact terms that contribute to 
this correlator. Converting the charge current given in ((3.19)) to momentum space 
we obtain 

ffli^P) = J (7^V4(^rn,.p' -p)YP-^E{u m ,p'). 

(3.24) 

After performing Wick contractions the correlator of interest is given by 

i(p) = —(j x (o,p)j y {o,-p)), 

IPz 

= +p) lf y S E (u m ,p)P + ]. 

(3.25) 
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Using standard methods to sum over the Matsubara frequencies we obtain 


E 


1 

(3 (iu m + ep~ 


tf(E p 
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1 

(3 ( iuj m + ep~ 


1 

tE p ) ) {ibJm, T T ^Ep-\-q) 

tep-) - uf(E p+q - uep _ - iu n ) 

“t“ tE p llEp+q 

iirfm + 6/1— 

6 / 1 — + i(jJ n uEp^-q^j 


Epf (Ep tc/l —) Ep+qf (Ep+q UCfl— 


E 


1 

/3 (iuj m + ep- 


tOJ n T tE p llEp+q 
(iu m + ep-) 2 _ 

tEp))(iu> m T ep— T l^n ^ Ep+q) 


tEpf(E p - ten-) - uE p+qf(E p+q - uep- - iu n ) 

tUJ n T tEp uEp_^q 


(3.26) 


While performing all these sums, we have ignored terms which are independent 
of temperature and chemical potential. Using the results for the sums as well as 
performing the trace in (3.25), we obtain 

e 2 f d 3 p' t tf(E p i +p - ten-) ~ uf(E p > - uep-) 

P ~ t ~± 2 J (2 n) 3 E p i +p tE p i +p — uEpi 

(3.27) 

After changing the variable of integration in the first term of the above equation, we 
obtain 

U ) = V — f d3p ' ( — E P'~ te 'tzl _i_ tuf(E p ,-uen-) v 

V ~ ( “ t 2 7 (2ir) 3 ^E p/ (tEpi - uEp, +p ) Ep> +p (uEp> - tE p , +p y 

(3.28) 


Summing over u in the first term and over t in the second term, we get, 


I(Pz) = 2e 2 


d 3 pi tf(E p ' - ten-) 
(2vr) 3 (E 2 - E 2 +p ) ’ 


2e 2 


(2U 2 t ± 


yy / dp'p' 2 sin 6tf(E p i — ten-) 

+ -L J 


(p 2 + 2p'p z cos 9) 


(3.29) 


At this stage we need to perform the angular integral. This integral can of course 
be easily performed, however we will show that we can take the p z —> 0 limit first 
and then perform the integral. Though the advantage of this procedure is limited 
in d = 4, we will see that for d = 6 this procedure simplifies the resulting integrals 
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considerably. Expanding in p z we obtain 
2e 2 


Hpz) = 


( 2 -) 2 ^ 


/ dp'p' 2 sin 9tf(E p i — tep) 

4- -L J 


v cos d 4p /2 cos 2 $ 
After a change of variables 2 = cos 9, the required integrals are of the form 


+ 0(p 2 z ) 


Ji = 


sin 9d6 
cos 6 


f 1 dz 

l-i z ’ 


J 3 = 


r sin Ode 

/ 0 cos 2 # 


f 1 dz 


(3.30) 


The integrals are all on the real line, we make these integrals well defined by shifting 
the integrand to z —> z + ie. Then we obtain 


J\ — 0 , J'2 — — 2 . 


(3.31) 


Such a prescription for doing the resulting integrals renders each term in the expan¬ 
sion in p z finite. In fact it is easy to check that performing the angular integral first 
and then expanding in p z , gives identical result to expanding in p z first and then 
doing the integrals term by term using the ie prescription. From these results for the 
integral, we obtain 

cl 1 ’ = lini — (j x (p)j y (-p)) = 

0 ip z 


Finally using ( 3.7) to obtain the transport coefficient we get 

Cr (4) ’* = Ci 4) = (3-33) 

Let us now relate this parity odd transport coefficient with the coefficient of the 
coefficient of the 17(1) anomaly. The anomalous conservation equation for the 17(1) 
current is given by, 

V„/ = [3 c A F^F pX + c m R^R e a J ■ (3-34) 

The gauge and mixed anomaly coefficients ca and c m for a single Weyl fermion are 
given by 12 


d Pitf(E pi -tep_), 

t=± d 

(3.32) 


CA 24ir 2 ’ ° m 192ir 2 ' 

Comparing the coefficient of the anomaly and the transport coefficient, we see that 
chiral fermions obey the relation 



12 For example, see (2.39) of [12]. 


Cl 1 ’ = -6c A /i_. 


(3.36) 













This relation obeys the replacement rule relating the gauge anomaly to the transport 
coefficient found in [12]. Alternatively, the coefficients c A and c m can be read 
out from the anomaly polynomial in four dimensions , 

V d = 3 +i(F, R ) = c a F A F A F + c m F A tr (R A R). (3.37) 

Evaluation of 


P'+P 



Figure 6: Diagram contributing to 


We first evaluate the relevant two point function in the Euclidean theory, this is 
given by 



lim t (j x (0,p z )T Ty (0, —p z ))- 

Pz^r 0 %p z 


(3.38) 


This correlator is defined by (2.15) and from the expansion of the action in (3.15), 
we see that it does not involve any contact terms. The expression for the current 
in momentum space is given by (3.19)) . Converting the stress tensor in (3.18)) to 
momentum space we obtain 




d 3 p' 

(2vr) 3 


+ P)i,Yw'y + + e/r_))P_^ E (p')- 


(3.39) 


After performing the Wick contractions we get 

—(/(0,pJT^(0,-p«))= (3.40) 

^Vz 

2/^7 S / 7^F3 Tv ^ XSe (p' + Pz){R T ipy + l y {iu m + ep-))S E (p)P+}- 

h'VZ m J V / 

For convenience we break up the terms resulting from the the Wick contractions into 
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two parts, 1 1 and J 2 , which are given by 


h{p z ) = 


d 3 p' 


I ^^[7 I 7Y/-P+]i>i(p , + P)t,«(p') u ^ 


iSfiPz ^ J (2vr) 
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4(27t) 3 J cos 2 9 
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A t (iu) m ,p' +p)A u (iuj m ,p'), 


xe 
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sind 


4(27t) 3 ,/ cos 2 6* 


d9d(f> / dpptfip — tep,-) + 0(p z ). 


(3.41) 


Here, in the second and the last line of the above equation, we have performed the 
Matsubara sums, expanded the integrand in the external momenta p z and finally per¬ 
formed the angular integrals using the ie prescription defined earlier. Note that the 
angular integrals are considerable simpler than the ones done first in [8] to evaluate 
this correlator. Using these results, we obtain 


cb = 


lim —{j x {0,p z )T Ty (0, -p z )) = lhn (h(p z ) + h{p z )), 

Pz^ o tp z Pz— >-0 


(3.42) 


le 


47T 2 


dp'p'(f(p' - epJ) + f(p' + ep-) , 


ie , on 7r 2 T 2 . 

-n(eV H-). 

8vr 2V 3 ; 


(3.43) 


We can now relate the above Euclidean correlator, to the desired transport co¬ 
efficient using (3.7). This results in the following expression 




vr^T 2 
3 ; 


(3.44) 


Finally let us relate the two terms, that occur in the transport coefficient <U 4 \ 
to the two anomaly coefficients ca and c m , defined in (3.34) and (3.35). First 
parametrize following [12] as 

(i 4) = -3c A p 2 - + c 4d T 2 . (3.45) 


Comparing with the results from Kubo formulae (3.44), we obtain 

e 3 _ e 

CA “ 247T 2 ’ C4d ~ _ 24’ 


(3.46) 
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From (3.34), (3.35) and (3.37), we see that value of ca obtained using the Kubo 
formula for <Cj 4) agree. This is a consistency check on our calculations. Finally we 
obtain the relation of the 

cm = -87r 2 c m . (3.47) 

This relation was obtained using consistency of the Euclidean vacuum in in [12]. 
This direct perturbative calculation provides a check on this argument. Note that the 
transport coefficient has been evaluated earlier by [8] , here we have demonstrated 
that the angular integrals simplify considerably on taking the zero momentum limit 
first. In fact doing this leads directly to the moments of the Fermi-Dirac distribution. 

Evaluation of A[, 4 




The last transport coefficient of interest in 4 dimensions is A 2 • The correspond¬ 
ing coefficient in the Euclidean theory is given by the following two point function 



lim - (T TX (0,p z )T Ty (0,p z )) 

Pz^r 0 ip z 


(3.48) 


From (2.16) and the perturbative expansion of the action in (3.15), we see that there 
is a contact term due to S ^ from the action. Therefore we obtain 

A< 4) = A + B, (3.49) 


A 

B 


~ Pz))e ' 


The expression for the stress tensors in momentum space can be read out from the 
following expressions 

T Ji (P) = J ~ p ^ Tip '* + + e/i_)]P_7Mp'), 

m ' ' 

T }i(-p) = yy 22 f + p)[Y¥y + + ep-)}P„y E (p')- 

(3.50) 
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This contribution from the contact term B in (3.49) was absent in the calculations 
of [10] and was considered in the recent work of [22], We will see that only after this 
contact term is considered, we will obtain the result consistent with the replacement 
rule. Performing the Wick contractions and the Matsubara sums, the term A reduces 
to 

(T TX (0,p z )T Ty (0, —p z )) 


A = lim 

Pz~>0 

1 


167r 2 


IPz 

((efi-f + ir 2 r 2 fi_) . 


(3.51) 


We now need to compute the contribution due to the contact term. This is given by 


B = -~ £ 1 


( 2lr ) 3 


' l P Trb'“7*7VP_] (\p' t „A,(iw m ,j>') 
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t=d= 
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d 3 p' 


tf (E p i - ten A), 


167T 2 

1 

487T 2 


(2vr) 3 ' 

J dp'p' 2 (f(E p i - e/i_) - f(E pi + e/i_)), 

(( e / i -) 3 + 7 T 2 T 2 / i _). 


Combining the results of term A and B, from (3.51) and (3.52), we get 

~^ 2 ] = y^((e^-) 3 + e7r 2 T 2 /i_)- 

We can now use the relation (3.7) to obtain the transport result of interest 

A^ 4) = -A (4) 


(3.52) 


(3.53) 


2 ; 

-1 


127T 2 


{{epA) 3 + eTT 2 T~p_). 


(3.54) 


Let us now verify that this result is consistent with the replacement rule. Fol¬ 
lowing [12], we parametrize transport coefficient as 

A^ 4) = 2 (~c A n 3 _ + c 4d T 2 p_). 


(3.55) 


Comparing this parametrization with (3.54), we obtain 

^3 


ca = 


24tt 2 ’ 


C4d 24' 


(3.56) 


From (3.35), we see that the value of ca is that of a single Weyl fermion. Furthermore 
c A d is identical to that obtained evaluating c{ 4) in (3.46) and it also satisfies the 
relation 


C/^d 87T c m . 


(3.57) 
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Dimension 

Correlator 

Value 

d=4 

ch 



BA) 

sMeV- + T) 


Af 

12 1 7r 2 ((e/t_) 3 + e7r 2 T 2 p_) 


Table 1: Transport coefficients In 4 Dimensions 


This relation was found by the argument involving the consistency of the Euclidean 
vacuum in [12]. It is important to note that the contribution of the contact term B 
was important to arrive at these results. 

The transport coefficients involving chiral fermions in d — 4 are summarized in 
the table (1). We note that in [19], the contribution of chiral gravitinos to these 
transport coefficients was discussed. However gravitinos cannot be charged in flat 
space without violating the gauge symmetry, -A- + d^e. Therefore they do 

not admit any U{ 1) charge current which implies there is no contribution to these 
transport coefficients from chiral gravitinos 13 . 

4 Anomalies and transport in d=6 

In this section we evaluate all the parity odd transport coefficients occurring in 
d = 6, at the leading order in derivatives, using the same methods developed in the 
previous section. This analysis will serve as a check of the computation in [22], of 
transport coefficients determined only by the charge current gauge anomaly. This 
was done using different approach than the one adopted in this paper. Our analysis 
will also use the ie prescription, developed in the previous section, to perform the 
angular integrals. Therefore the results of the d = 6 calculations also test this 
prescription. Finally we also determine the contributions to the transport coefficients 
from the mixed as well as the pure gravitational anomaly, which was not done in 
[22] . Our results are in agreement with the replacement rule for the mixed and pure 
gravitational anomaly, found by [7, 12] for chiral fermions. But for the case of chiral 
gravitinos which contribute to the pure gravitational anomaly, we see that the their 
contribution to the transport coefficients does not obey the replacement rule. We 
precisely identify the terms which contribute to the anomaly polynomial but not to 
the transport coefficient. 

We begin with parametrization of the constitutive relations in d = 6, where we 
focus on the the leading parity odd transport coefficients. Following the parametriza- 

13 We thank Zohar Komargodski for emphasizing this to us. 
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tion given in [ 22 ], the charge current and the stress tensor is given by 

A6) 

nv? + ^-e^ s u u F aP F lS 

AG) AG) 

+ ^-e^ s UudaUf}Fl5 + ^e^ 6 u u d a u^us, 

\(G) 

(e + p)u»u v - pg^ + u^^-e^ s u p F^F j5 

A (6) a (6) 

+ AL_ e »P«W UpdaU p FjS + FL_ e ^s UpdaU ^ Us) + ^ 

(4.1) 

where, e,p,n,p are the thermodynamic variables denoting energy density, pres¬ 
sure and charge density. The six anomalous transport coefficients are denoted by 
A^, Ci 6 \ C 3 6 " 1 • The Minkowski metric, we work in, has the signature 
(rj /iu ) = diag (1, —1, —1, —1, —1, —1). Let us also label the directions as follows, 
the position vector or momentum vector has components in directions labeled by 
(t,a, z,x,y,b). We choose the orientation in d = 6 so that e tazxyb = 1. To arrive at 
the Kubo formula, we assume that the fluid is perturbed from rest and its velocity 
field is given by u y = (1,0, v z : v x , 0, 0)) . The velocity field v z (t,b) depends only on 
time t and coordinate b, while the velocity field v x (t,z ) depends only on time and 
the coordinate z. We then consider the metric perturbation h tz (t,b ) and h tx (t,y), 
while the rest of the fluctuations are set to zero. Similarly we consider gauge fluc¬ 
tuations a z (t,b) and a x (t,y). We summarize the perturbations and the coordinates 
they depend on in the following equations 

u y = (1, 0, v z (t,b), v x (t,z), 0, 0), (4.2) 

u p = ( 1 , 0 , - v z (t , b) + h tz (t , 6), -v x (t, z) + h tx {t , z), 0 , 0 ), 

A p = (/i, 0 , a z (t,b), a x (t,y ), 0 , 0 ), 

g»v = Vfiu + Vi with h tz{t, b ), htx(t, y) ^ 0 . 

With this choice of the polarizations and perturbations, it can be easily seen that the 
contributions proportional to the anomalous transport coefficients in (4.1), appear at 
second order in fluctuations and derivatives. Therefore the Kubo formulae for these 
transport coefficients are given by three point functions. Expanding the constitutive 
relations to second order, we obtain following expressions relating the three point 


f = 


rpilV _ 


40 






(4.3) 


functions to the transport coefficients 


8 ^( 6 ) 

5v(6) 

3^,2 

2 >( 6 ) 

3 ^ 3 


8 

3 

2 

3 

2 

3 


A< 6) 

A? } 

Af 


= — lim 

ky,p b ^0 


= — lim 

ky,Pb^0 


= — lim 

ky,Pb^> 0 


= — lim 

k y ,P6-A0 


= — lim 

ky,Pb ~>0 


= — lim 

ky,Pb^ 0 


(j a Q + k)j x \-k)j z (-p)) 
ikyipb 

(j a (p + k)j x (—k)T tz (—p )) 

(j a (/c + p)T tx {—k)T tz (—p)) 
ikyipb 

( j a (p + k)j x (—k)T tz (—p )) 

( j a (p + k)T tx (—k)T tz (—p)) 

ikyip b 

(T ta (p + k)T tx (—k)T tz (—p)) 
ikyipb 


To un-clutter our notations, it is understood that the external frequencies are set to 
zero from the beginning and that the only non-zero component of the momentum 
vectors p,k are in the b and the y directions respectively. From (4.3), we arrive at 
the following relations between the transport coefficients 


4Af> = Cf, Af = Cf. (4.4) 

Therefore it is sufficient to focus on evaluating the coefficients ( 3 ^ and \ f ^. 

The strategy for evaluating the correlators will be same, we will first evaluate them in 
the Euclidean theory and then analytically continue them to the Lorentzian theory. 
Let the value of the corresponding correlators in the Euclidean theory by , ( 3 ^ 

and A^. Then going carrying out the analytic continuation leads to the following 
relations 


c?> = cf, cf = <cf , ( 4 . 5 ) 

cf = -cf. Af = -nf ■ 

We will now proceed to evaluation of the contribution of these correlators, first 
for chiral fermions and then the contribution of the chiral gravitinos to A 3 6 ). 

4.1 Chiral Fermions 

For Weyl fermions in d = 6 , the partition function is given by 

j V^Viiexp(S E ), (4.6) 

where ^ is a Euclidean Dirac spinor in d — 6 and the Euclidean action is given by 

S E = J dTd 5 x y /ge^'fD li P-.rj>. (4.7) 
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The covariant derivative is defined as 


D p i) = dp'i/j + i u pab a ab ij + ieApijj, (4.8) 

where, uj pa j, is the spin connection and a ab = The gamma matrices follow 

the same conventions as before. The conventions for the metric, vierbein and their 
perturbations remain same as that in d — 4 which is given in equations (3.11), (3.13) 
and (3.14). The chiral projection operator is defined by 

P± = ^(1 ± 7 ? )> 7 7 = iYl a Yl x l y l b - (4-9) 

To evaluate the Feynman diagrams, we will need the following trace 

Tr(7 P 7V7V 7 V) = -8e pqrstu , e Tazxyb = i. (4.10) 

Wick contractions are performed by the propagator given in (3.21). Expanding the 
action in in these fluctuations, we obtain 

= J drd 5 x^(YDr + 'y i d i )P-'ip, 

= - J drd b x h d^^{Yd v + Yd p )P_^ 

+ie J drd 5 x^ r f ‘+ 0(ha ) + 0(/? 2 , a) + 0(h 3 ), 

= [ drd 5 x hxadflhua ^T pXl/ P_rj) + 0(h 3 ). (4.11) 

J 16 

At first glance one would have expected that we must expand the action to at 
least 0(h 3 ),0(h 2 a ) to obtain all contributions to the correlators of interest. This 
is because now we are evaluating three point functions. However terms involving 
0(h, a ), 0(h 3 ), 0{h 2 , a) in the expansion of (4.11), do not have sufficient num¬ 
ber of T matrices to ensure a non-zero result for the correlators of interest. We 
need 6 T matrices to ensure non-zero result in d = 6. There is another reason the 
0(hA) terms do not contribute. From the structure of the perturbations in (4.2) 
of the metric and the gauge field, it is easy to see there are no terms of the type 
0(h TX a z ),0(h TZ a x ) and therefore there is no contributions to the correlators of in¬ 
terest. Finally the reason 0(h 3 ) term is not considered in S^ of (4.11) is again 
due to the fact that such terms do not have sufficient number of gamma matrices 
to saturate the correlator of interest. Therefore the terms in the expansion (4.11) is 
sufficient for evaluating all the transport coefficients of interest. 
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Evaluation of 

We begin by evaluating the coefficient C[ 6 \ the corresponding correlator is given by 

1 


O'W> = 


z 


(o) 


5 2 A 


v^v^' ( SSe 6 -A SSe 

V V 9 $ A V V 9 $ a „ V 9 $A p 

S 2 S e 


(4.12) 


5S f 


+ 


5S f 


y/g^Apy/g8A v y/gbAp y/gSA^SAp y/gSA, 

5S e S 2 S e S 3 S e 


y/gSAvy/gSAvy/gSAp, 

Here it is understood that only the connected correlators should be evaluated while 
performing Wick contractions. The correlator of interest is given by 


C (6) = - lirn 3 ti a ( k + P)j x {-k)j z {~p)) 

1 Pb,ky ~>0 8 


(4.13) 


ipbiky 

As we have mentioned earlier, we choose the kinematic configurations in which the 
momentum vectors p, k have components only in the b, y directions. The frequency 
dependence of the external momenta has been set to zero from the start. First 
note that from the expansion of the action in (4.11), the gauge fluctuations appear 
linearly and therefore there are no contributions from terms involving 2 derivatives 
of the action in (4.12). Therefore the correlator in (4.13) is obtained by performing 
Wick contractions of the flat space charge current which is given by 

d 5 p' 




(2tt) ! 


:^\uj m ,p' - p)YP~^{^ m ,P')- 


(4.14) 


After performing the two possible Wick contractions, we obtain 

(j‘(p+k)f(-k)f(-p)) = (4.15) 

m ' ' 

x (Tr^S (u} m , p + k + p)^/ x S (u> m , p + p)YS(u m , p)P+] 
+Tr[y z S(u) m ,p' - p)^ x S{uj m ,p' -p - k)-f a S E (um,p)P+}) ■ 


The two Wick contractions are summarized in the figure 8. The details of carrying 
out the Matsubara sums and then the angular integrals using the ie prescription is 
given in the appendix B.l. The end result for the Euclidean coefficient is 

o p 4 roc 

Cl 6) = {f (p - e/x-) - f{p + e/i_)), (4.16) 

3e 4 

647r 3 ^ 


Using the relations for analytic continuation given in (4.5), we obtain 

3e 4 


c! 6) = 


647 T 3 


p— 


(4.17) 
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Figure 8 : Diagrams contributing to 

Let us now compare these coefficients to that in the anomaly polynomial in d = 6. 
The anomaly polynomial for a single Weyl fermions in d — 6 is given by [19] 

v d . w (F, R) = cj‘ + cpF 2 T?(R 2 ) + c,(Tr (Ft 2 )) 2 + cj^Tr/fl 1 ) - gTrtfl 2 ) 2 ), 

(4.18) 

where 


C(y. 


192tt 3 ’ 
7 


c 7 = 


1843207T 3 ’ 

Thus the transport coefficient satisfies the relation 

Cl (6) = -9 Call— 


C/3 - - 

7687T 3 ’ 

1 


11520tt 3 


(4.19) 


(4.20) 


This relation is the expected relation between the gauge anomaly coefficient and the 
transport coefficient of interest, predicted by general considerations of hydrodynamics 
in [5, 12]. 

Evaluation of 

The correlator of interest for the transport coefficient is given by 


( ffTO*) = 


'D^T>Tpe SE x 


Z 


(o) 


(- 


5S r 5Sp 


5S f 


+ 


6 2 S f 


5Sf 


\Sy/gA p ^g5A v ^fg5g p<J y/gSA^y/gSA,, y/gSg^ 
| 8 2 S e 5S e | 5S e 5 2 S e 
y/gdA^yfgdgpe JgdA v ^g5A p ^g5A u 5^gg pa 
S 3 S e 


y/gSapy/gSAv^/gSg pa 


(4.21) 
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Here again it is understood that only connected Wick contractions need to be eval¬ 
uated. After evaluating the correlator in Euclidean space, we extract the following 
coefficient 


y(6) _ _ 

S2 ~ 


Hm 3 ( j a (k + p)j x (-k)T TZ (-p )) 
m \o 2 


(4.22) 


Pbiky —>0 2 ZPb'lky 

From the expansion of the action in (4.11), we see that none of the contact terms 
contribute to this correlator and therefore we can evaluate the correlator in (4.22), 
by performing Wick contractions on the flat space charge current given in (4.14) as 
well as the stress tensor given by 


T ji z (~P) = f + P)^ Ti P'z + r( iuj m + 

(4.23) 


The two possible Wick contractions are summarized in the figure 9. The details 




Figure 9: Contributions to 


of the Wick contractions, Matsubara sums and angular integrals resulting in the 
calculation are provided in B.2. The end result of the corresponding Euclidean 
transport coefficient is given by 


y(6) 

S2 


3e 2 

i 16T^J p ^^ p ~ ep ~"> ~f(P + e P-)i 


. 3e 2 
1 32n 3 


(eV_ + 


3 ’ 


(4.24) 


Again using the relations in (4.5), we obtain the transport coefficient of interest to 
be 


y(6) 

S2 


3e 2 
32t r 3 


(e 2 /i_ + 


tt 2 ^ 
3 ; 


(4.25) 


Let us compare this coefficient with the coefficients in the anomaly polynomial. 
For this, following [12], we first parametrize the transport coefficient as 


A 6 ) 

S2 


3(—6 Cq/i 3 . + c^T 2 ), 


(4.26) 
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where c a for Weyl fermions is read out from (4.18). Comparing (4.25), we see that 
c .^ is related to the mixed gravitational anomaly coefficient cp by 

c“ = -8n 2 c p . (4.27) 

In [12], such a relation was obtained using the argument involving the consistency 
of the Euclidean vacuum. Here we have explicitly demonstrated this relation using 
the perturbative computation. 


Evaluation of 


Now consider the correlator 




(fr f T n = -W / 

'D^'D^e s<E ' 1 x 



( SSe 

SS E 5S e 

5 2 S e 

SS e 


Ky/gSAp ^fgSg v p y/gfigpa y/pSA^y/pSgvp y/dSgpa 

5 2 S e SS e 5S e S 2 S e 

yfpSAp\/f]3fjprj A J~pSg v p y/HSAp \fQ^9vf3\[Q^Qpa 

+__y 

y/goAfty/gdgupy/gdg pa J 

The corresponding transport coefficient in the Euclidean theory is given by 


7(6) 

83 


— lim 

Pb,ky~>0 


3 (j a (k + p)T TX {-k)T TZ (-p )) 

2 ipb'iky 


(4.28) 


(4.29) 


Now from examining the expansion of the action in (4.11), we see that the only 
non-zero contributions to the correlator are given by 


(f T "P T P°) = =A + B, 
4 


(4.30) 


A = 


B = 


Z 


(o) 


'Difj^T>if>e kJE 


(o) SSp 


SSf 


5S f 


( 0 ) 


z 


(o) 




y/gSAp y/g5g uP y/gSg pa 

ss E 


5 2 Sf 


y/gSAy y/gSgvpy/gSg pa 

Therefore we have a contribution from the Wick contractions of the flat space currents 
denoted by A, together with a contact term which involves Wick contraction with a 
flat space current denoted by B. These are summarized in figure 10. and figure 11. 
respectively. We have evaluated these contributions in detail in appendix B.3. The 
end result for the Euclidean correlator is given by 


7(6) 

83 


-3e 


167T 3 
—e 
167T 3 


dpp 2 (f(pi - e/i_) - f(p + eg ,-)), 


[(e/i_) 3 + eix 2 T~g_\. 


(4.31) 
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Figure 10: Contributions to 


from Wick contractions of flat currents. 


p'+p+k 



h 3 h h 

IX b TZ 



Figure 11: Contributions to Q 6 


from contact terms. 


Again using the relations in (4.5) to go from the Euclidean correlator to the transport 
coefficient of interest in the Lorentzian theory, we obtain 

<f ,H = y^K^-) 3 + ^ 2 TV-]. (4.32) 

Let us compare this transport coefficient to what is expected from the replace¬ 
ment rule. Following [12], we parametrize the coefficient as 

$>•* = 6(-2c«/i 3 _ + 5“T 2 /i_). (4.33) 

Comparing the final result of the perturbative evaluation in (4.32) and the above 
equation we see that indeed the results for the coefficient c a agrees with the gauge 
anomaly of a single Weyl fermion given in (4.18). Also we have the relation 

= -8vr %, (4.34) 


which is again consistent with that obtained by the argument involving the consis¬ 
tency of the Euclidean vacuum. 





Evaluation of 

The last transport coefficient in d — 6 is given by the following correlator 


^ rj-lHCtrj-aj(3rppcr ^ _ 


5?(0) / 

Z, E J 

( SS E 


V^Vipexp x 


(4.35) 


5S f 


5S f 


+ 


S 2 S f 


SSf 


\yfgSg f.ia y/gfigvp VgSg pa yfgSg flOL y/gSg„p a /gS pa 

| S 2 S E SS e SS e S 2 S e 

\fg^g pa \fg^g pa \fgSg v ^ yfgSgpa \fg^gvfi\fg^ pa 

S 3 S e 


yfgSg fia y/gSgufiy/gSg pa 

The coefficient in the Euclidean theory is given by 


= — lim 

2 Pb,ky^0 


-(T ra (k + p)T TX (-k)T TZ (-p)) 

ipbiky 


(4.36) 


From the expansion of the action in (4.11 ), we see that only non-zero contributions 
to this correlator arise from 


(T^T^T^) = A + B, 
8 


A = 


B = 


Z 


(o) 


Z 


(o) 


5 2 S. 


V^V^exp (S { E ] ) 

V^Vipexp {S { E ] ) 
SS E 


SSf 


SSf 


5S f 


VgSg f.ia VgSg»e \fgSg pa 

e SSe 


S 2 Sf 


+ 


yfgSg /ia VgSgvp y/gS pa 
SSe S 2 Se 


(4.37) 


y/gSgpa^/gSg p(7 yfgSg u p y/gSg^a \fgSg v ^^jgS t 


pa 

Term A involves Wick contractions of the flat space stress tensor while B is a contact 
term. From the explicit evaluation of term B in appendix B.4, it is seen that all terms 
in B contribute equally. 

The contributions to A and one of the terms in B are summarized in figure 12 
and 13 respectively. The detailed evaluation of all the diagrams yields the following 
result 


\( 6 ) _ 

As “ 16vr 3 


-3 i 


dppl(f{E p i - e/i_) - f(E pi + e/x_)), 


(4.38) 


-3 i , 7n 4 T 4 


i o Z ZrriZ Z , 4 4 \ 

+ 2e 7 r T fi_ + e fi_) 


647r 3 15 

Writing down the transport coefficient in the Lorentzian theory using the relations 
in (4.5) gives 


\(6) _ 3 7n 4 T 4 # n ^2^2rr2..2 

3 64tt 3 ^ 15 


+ 2e z 7T z T z pi+e 4 pi). 


(4.39) 
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Figure 12: Contributions to 


from Wick contractions of the stress tensor 




Figure 13: 


Contributions to A 


( 6 ) 

3 


from contact terms 


The importance of this transport coefficient is that it captures one of the pure grav¬ 
itational anomaly d = 6. Again following [12], we parametrize this coefficient as 

a£ 6) = 9 (—c q /t 4 _ + c“ry_ + cf T 4 ). (4.40) 


Comparing the result of the perturbative evaluation in (4.39), we see that the values 
of the gauge anomaly c a and the coefficient c^ is related to the mixed anomaly, as 
seen in the earlier correlators, while c 6 g d is given by 


g6d_ 77r 
9 2880’ 


(4.41) 


Comparing with the coefficients of the pure gravitational anomaly of a Weyl fermion 
c 7 given in (4.18), we obtain 


c 6d = -(8 t r 2 ) 2 c. 


(4.42) 


This relation is what is expected from the replacement rule and the consistency of 
the Euclidean vacuum 14 . 


4.2 Chiral Gravitinos 

Chiral gravitinos cannot be consistently coupled to a charge current. Therefore 
they contribute only to the pure gravitational anomaly and only to the transport 
coefficient X-f' 1 . In this section, we show that analogous to the coefficient A 1 - 2 ) for 

14 We have not checked if the precise coefficient of — (87t 2 ) 2 agrees with the normalizations of [12]. 
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the chiral gravitino like theory d — 2 in section 2.3 , the for chiral gravitinos in 
d = 6 do not obey the relation c 6 g d = — (87t 2 ) 2 c 7 . We will identify the terms which do 
not contribute to the transport coefficient but contribute to the anomaly polynomial. 
The gauge fixed gravitino Lagrangian is given by 

S E = S£ + Sf, (4.43) 

s e = I dTdx 5 y /ge£il>lD li P-.il) p , 

where the covariant derivative is given by 

D,r = d^ p + \oj^a cd r + r (4.44) 

The action S^ 1 is the action of the corresponding ghosts. The contribution of all the 
ghosts can be accounted by subtracting the contribution of a chiral spinor. The action 
S£ is similar to the action of 6 chiral fermions except that the covariant derivative 
in (4.44) contains the Christoffel symbol in addition to the spin connection. The 
linearization of the Christoffel symbol term results in the following interaction term 
in addition to the terms of the type present in (4.11). 

^Christoffel = \ J dTd b x y /ge^ p 'frf a {d lx h aa + d a h aix - dji^ + 0{h 2 ))^ c . 

(4.45) 

The linear terms in the metric modify the stress tensor. There are higher order terms 
in the expansion of the metric and one would have expected possible contributions 
to contact terms from them. However such terms do not have sufficient number of 
gamma matrices to contribute to the correlator involving 3 stress tensors. We need 
6 gamma matrices, so that the correlator is non-zero. Any possible contact term 
from higher order terms in the above action have lesser number of gamma matrices. 
Therefore we conclude that the only difference of the chiral gravitino, compared to 
that of the chiral fermion, is that its flat space tensor acquires additional term, due 
to the linearization of the Christoffel symbol. The stress tensor is given by 

T pv = Tg + T^, (4.46) 

Tg 2 = -^(9 (r ('0 t ’ CT (x)7 - d a (4} p \x)YP-V{x) + (/i ^ v)). (4.47) 
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Note that the contribution Tjj% is due to the linearization of the Christoffel symbol. 
In momentum space the stress tensor is given by 


Tjm = Tlf(k) + T™(k), 
d 5 p' 


(4.48) 


Tgi(k) = I - fc )(7 T iPa + 7 x iUm)P-V(u m ,p'), 


Tg2(k ) = - Ph T P-^ S (Um,p) - ^ T \u mi p - p)T P-V (u mi p') 


2/3 ^ (27r) 5 

m v 7 

+ (r V-)- s)). 

We can now evaluate the correlator 


Ag 6 - = — lim 

2 Pb,ky^> 0 


■(T Ta (k + p)T TX (-k)T TZ (-p)) 

ipbiky 


(4.49) 


Apart from the contact terms arising from the linearization of the spin connection, 
we have to evaluate the following contribution to the correlator 


A=(T?{k + p)T?(-k)T5r(-p)). 


(4.50) 


Substituting the expression for T// s from (4.48), it can seen from the structure of the 
Wick contractions, that the only non-zero terms are 


A — A\ + A 2 (4.51) 

M = (T T g ?(k+p)T™(-k)T™(-p)) 

A 2 = (T£(k + p)T™(-k)T™(-p)) + (T; 2 a (k + p)T™(-k)T£(-p)) 

+ (T£(k + p)T™(-k)T™(-p)). 


After performing Wick contractions, the sum of the terms in A 2 organize themselves 
to 


Ao — —— 


d 5 p' 


4/3 ^ J (27t) 5 (4 ' 52) 

(p^Tr [y a S(uj m: p' +p +k)(YPx + l x iu m )S{u m ,p' + p)YS(u m ,p')P + ] 

+Pb Tr ['fS(u rn ,p l - p)(Yp' x + l x iUm)S(u m ,p' -p - k)-f a S(u m ,p)P + } 
+kyTr [^ a S{u m ,p' +p+ k)^ x S(u m ,p' +p)(Yp' z + 7 z iUm)S(u m ,p , )P+\ 

+k 2 y Tr [(Yp' z + l z ^m)S{u m ,p' - p)'y x S(u rn ,p' - p - k)^ a S{u m ,p')P + \) . 


After inserting the propagator and performing the Matsubara sums, it can be seen 
that each of the terms in A 2 is proportional to pbk y . Therefore, the leading contri¬ 
bution to the transport coefficient from A 2 is of the form 


Ao 


lim 

Pb)ky —^0 ZpbZrvy 


— 0(pl) + 0{k 2 ) 


(4.53) 
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The entire contribution from A 2 to the transport coefficient vanishes in the zero 
momentum limit. This is the same mechanism by which, the contribution from terms 
arising from the linearization of the Christoffel symbol to A*- 2 ' in the gravitino like 
theory in d — 2, vanishes. The contribution from A\ is 6 times that of a chiral fermion. 
The contribution from the contact terms arising from the linearization of the spin 
connection is also 6 times that of a chiral fermion. To these contributions we must 
subtract the contribution of a single chiral fermion to account for the ghosts. Using 
the results of the the section 4.1 equation (4.39), we conclude that the contribution 
of a chiral gravitino to A 3 is given by 


A (6) 

gravitinos 


3 7vr 4 T 4 
647r 3 15 


qp 6d 

gravitinos* 


(4.54) 


Here the last equality arises from the parametrization following from the replacement 
rule. 

Now let us compare the transport coefficient to the coefficient of the pure gravi¬ 
tational anomaly. The terms arising from the linearization of the Christoffel symbol 
do give a contribution to the anomaly polynomial of the gravitino [23] 15 . In fact 
such terms give rise to the interaction of the magnetic moment of the gravitino with 
the gravitational field. These terms modify the gravitational anomaly coefficient c 7 
for gravitinos which is given by 


C7,gravitinos 


-275 x 


1 

184320vr 3 ' 


(4.55) 


Comparing (4.54) and (4.55), we see that gravitinos do not obey the relation c® d = 
— ( 87 t 2 ) 2 c 7 , which is argued from the consistency of the Euclidean vacuum in [12]. 
Here we have identified the terms which contribute to the anomaly polynomial but 
not to the transport coefficient. From our analysis in d — 2, d = 6 we see that 
gravitinos contribute to correlators involving only the stress tensor, since they are 
not charged. Their contribution is essentially d — 1 times that of a chiral fermion to 
the transport coefficient. This is the same conclusion arrived at, by analysis of the 
hclicity of the thermal state, in [19]. 


15 See equation (35) of [23] 






Dimension 

Correlator 

Species 

Value 

d=6 

cr 

Fermions 



716 ) 

S2 

Fermions 

Jg,(eW + +7 


7 ( 6 ) 

^3 

Fermions 

1^3 [(e/i-) 3 + e7T 2 T 2 /i_] 


\T61 

A 3 

Fermions 

s y(+p+2 e w( M _) 2 + e V-) 


\(S) 

A 3 

Gravitinos 

15 (7tt 4 T 4 \ 

64-7T 3 ^ 15 ' 


Table 2: Transport coefficients In d=6 


5 Hydrodynamic modes in d = 2 Weyl gas 

A charged 1 + 1 dimensional relativistic fluid admits 2 hydrodynamic modes, the 
sound mode and the charge diffusion mode. The leading dispersion relation for these 
modes, without anomalous terms in the constitutive relations for the stress tensor 
and the charge current, was obtained in [24], In this section we would like to show 
that due to the presence of anomalous terms, proportional to and A^ in the 
constitutive relation, the charge diffusion mode acquires a non-zero speed. We then 

apply the discussion to the case of ideal d = 2 Weyl gas, held at a chiral chemical 

potential p_. As seen in section 2 , such a gas has a non-zero value for and A^. 
The speed of sound for this gas remains the same as the speed of light, however the 
speed of the charge diffusion mode is non-zero and is equal half the speed of light. 

Consider a relativistic fluid in d — 2 with anomalous transport and A 1 - 2 ) 
non-zero. The constitutive relations to the first derivative in velocity in the Lorentz 
frame is given by 

= (e + P) v?u v - Prf v - £{u»u v - rf v )d x u x , 
f = pu » - aT{rT - u*u v )d v ^) + C ,^u v . (5.1) 

The signature of the metric is given by rj lw = diag(l, —1). u M is the velocity field 

with — 1 , £ is the bulk viscosity, a the conductivity and Q is the anomalous 
transport coefficient, e, P, p, fi refer to the energy density, pressure, charge density 
and chemical potential respectively. Q. in the Lorentz frame is related to the transport 
coefficients in the anomaly frame by 

<« = C (2) - (5.2) 

To study linearized hydrodynamic modes, consider small fluctuations from the rest 
frame with velocity field (l,v 2 ) ( \v z \ « 1). We also have fluctuations in the energy 
density, pressure which are related to fluctuations in the stress tensor and charge 
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current, using the constitutive relations in (5.1). These are given by 


tf) xpOz 

ST 00 = Se, 8T 0z = (e + P)8u z , ST ZZ = SP - * , 

v i e + p 

8j° = 5p — ( v 8u z , Sj z = p8u z — aTd z (Sp) — 5( v . (5.3) 

It is convenient to use the energy density e and charge density p as independent 
variables. We can write the variation in chemical potential and pressure as 

dp = d e p5T 00 + d p p5p, SP = d e PST 00 + d p P8p. 

where p = p/T. The hydrodynamic modes are eigen modes of the linearized equa¬ 
tions of hydrodynamics, which are given by 

d tJ ST' w = 0, d p 5f = 0. (5.4) 

Here we set the background gauge field strength ,including its fluctuations, to zero. 
Writing out these equations in Fourier space 16 , eliminating 5T ZZ , and 5u z using (5.3), 
we obtain the following set of equations from the conservation equation for the stress 
tensor 


- iuST 00 + ik5T 0z = 0, (5.5) 

ST 0z 

- iuST 0z + ik(d e PST 00 + dpPSp) + = 0. 

Now using the first equation in (5.5), to eliminate ST 0z , we obtain 

Sp(dpPik) + ST 00 (—^ + ikd t p] = 0. (5.6) 

\ k e T P J 

A similar analysis of the conservation equation for the current yields 


Sp(—ico + k 2 aTd p p — ikdp( v ) 
+5T 00 f k 2 aTd e p + 


(5.7) 


00 I *-™r. i iU P - ikd'Q- ^ 


e + P 


k(e + P) 


= 0 . 


From the homogeneous equations (5.6) and (5.7), we see that to obtain a non-trivial 
solution for the fluctuations Sp, ST 00 we must have the condition 


(' dpPik ) ( k 2 aTd e p H——^ — ikd e £ v + ^ 


ojk£ 


e + P ' k(e + P) 

+ ikd e PJ {—iu + k 2 aTd p p — ikd p ( v ) 

16 For eg. Fourier mode of ST ZZ = e~ iut+ikz 6T zz (u,k). 


k 


+ 


e + P 


= 0 . 


(5.8) 
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To obtain solutions which satisfy the condition in (5.8), we assume the dispersion 
relation of the form 


u — vk — iDk 2 + 0(k 3 ). (5.9) 

Substituting this anstaz into (5.8), we see that the velocity v is determined by the 
following cubic equation 

v 3 + v 2 (y~ d ^p + d p ( v ^J - v (^~j 5 + deP^J ~ d e Pd p ( v + d p Pd e ( v = 0. (5.10) 

This equation is obtained by setting the coefficient of the 0(k 2 ) term to zero on 
substituting the anstaz (5.9) in (5.8). As a simple check of this equation, consider 
setting the anomalous transport coefficients ( v = 0 , then the equation reduces to 

v 1 — + deP^j v = 0. (5-11) 

The three solutions to this equation are 

V = ± (f^ + ^ P ) 2 ’ v = 0 - (5-12) 

The first two roots correspond to the speed of sound and the last one is the speed 
of the charge diffusion mode. Now one can show for ( v ^ 0 but small, there are 
corrections to both speed of sound and the speed of the charge diffusion mode that 
depend on ( v . On further analysis, one can show that the the dissipative constant D 
for both the sound and the charge diffusion mode receives corrections on turning on 

Cv 

We will now show for the ideal Weyl gas the three roots of (5.10) are v = ±1, v = 
—1/2. The thermodynamics of the ideal Weyl gas is given by [18] 


e 

P 


e 2 /i_ T 2 7 T 
4tt + ~12~" 
e 2 /i_ 

2vr 


P, 


(5.13) 


Using (5.2) and the values of and we obtain the transport coefficient to be 
given by 


Ct, = 


e 2 /m 
2t r 
P 
2 ' 


P 


Air P + e 


7T 2 r 2 


2 2 i 

e [i_ + 


(5.14) 


To obtain the last equality we have used the expression for energy density and e = P. 
Since the ideal Weyl gas is conformal and there are no interactions the bulk viscosity 
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£ as well as the conductivity a vanishes. Now from the relations in (5.13) and (5.14) 
we see that 


d p P = 0, d e P = 1, 



(5.15) 


Substituting these values in equation (5.10) for the velocity of the modes we obtain 


, 1 9 1 

v + -v -v -= 0. 

2 2 


(5.16) 


Therefore the 3 roots are 


1 


(5.17) 


- = ± 1 . ~2 


As we discussed earlier, the first two roots correspond to the sound mode which 
moves with the velocity of light. The last root corresponds to the speed of the 
charge diffusion mode. Note that there is a unique sign for this speed, this is because 
of the definite chirality of the Weyl fermions. If the we had considered Weyl fermions 
with opposite chirality the sign would have been positive. 

It will be interesting to explore this phenomenon further to see if it is observ¬ 
able in realistic low dimensional quantum systems like quantum wires. It will be 
also interesting to study the effect of the anomalous transport coefficients for other 
relativistic systems in 1 + 1 dimensions. One such example is the charged Dl-brane 
which was studied in [24]. 

6 Conclusions 

We have evaluated all the leading parity odd transport coefficients for chiral fermions 
in d = 2,4,6 as well as chiral bosons in d — 2 at one loop. The analysis was done 
using finite temperature held theory. We have seen that chiral fermions as well as 
chiral bosons in d — 2 obey the replacement rule of [12] which relates the coefficient 
of anomalies including both mixed and pure gravitational anomalies to the respective 
transport coefficient. As far as we are aware the calculations in d — 2 is the first 
direct perturbative evaluation of transport coefficients using the Kubo formula. The 
calculations in d = 4,6 serve as an alternative approach as well as a check of the 
calculations done in [22] which involved Schwinger-Keldysh propagators. We have 
also kept track of the contributions from the gravitational anomalies which was ig¬ 
nored in [22], For a chiral gravitino like theory in d — 2 as well as chiral gravitinos 
in d — 6 the relationship between the pure gravitational anomaly and the corre¬ 
sponding transport coefficient found using the argument involving the consistency of 
the Euclidean vacuum breaks down. We have identified the mechanism of how this 
occurs in the perturbative analysis. We also developed a simple method to perform 
the angular integrals that occur in the transport coefficients which involves taking 
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the zero momentum limit first and the use of an it prescription. Finally we have 
shown that the charge diffusion mode of the ideal Weyl gas in d — 2 at finite chemi¬ 
cal potential acquires a speed which is 1/2 the speed of light due to the presence of 
anomalous transport coefficients. 

There are several directions in which the work done in this paper can be extended. 
It will be interesting to extend the analysis done in this paper to all even dimensions. 
A more challenging task would be to evaluate the contribution of self dual tensors 
in d — Ak + 2 for k > 0 to the transport coefficient Though the coupling of these 
tensors to the charge current is not explicitly known, one can evaluate contributions 
related to the pure gravitational anomaly of these tensors. Direct holographic checks 
of these transport coefficients in AdS 3 , AdSi backgrounds similar to the ones done 
in AdS$ in [9] is another direction direction of interest. Finally it will be interesting 
to study to what extent these parity odd transport coefficients are one loop exact 
analogous to what was done for the coefficient in d — 4 by [11]. 
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A Moments of statistical distributions 

To evaluate the various transport coefficients of interest in this paper we require 
several moments of the Fermi-Dirac distribution. Let f(p) be the Fermi-Dirac dis¬ 
tribution which is given by 



(A.l) 


The results for the following moments are used at several instances in the paper 



o 


■OO 


dpp 3 [f{p 


e/i) + f(p + e/z)] = + 2 eVTV + e 4 //) (A.5) 
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These integrals are derived in [18]. Let b(p) by the Bose-Einstein distribution which 
is given by 

b (P) = ( A - 6 ) 

Then the following integral is required in the analysis of chiral bosons. 


dpp b(p) = 


7 r 2 T 2 


(A.7) 


B Evaluation of transport coefficients in d = 6 

In this appendix we detail the calculations involved in obtaining the leading parity 
odd transport coefficients in d = 6. 

B.l Contributions to 

Before we begin, let us recall the kinematic set up of the correlator considered in 
(4.13). This kinematic configuration will be the same for all the other correlators 
considered in d = 6. There are two external momenta p, k. We set the frequency 
components of both these external momenta to zero. The momentum vector p has a 
non-zero component only along the b direction while momentum vector k has com¬ 
ponent only along the y direction. While performing the integration over the spatial 
directions of the internal momentum labelled as p' it is convenient to parametrize its 
components in terms of angular variables as 

p' b = |p'I cos 01 , p'y = |p'I sin 01 cos 02 , Px = |p'I sin 01 sin 02 cos 03, (B.l) 

p z = |p'| sin 0 i sin 0 2 sin 0 3 cos 04 , p' a = |p'| sin 0 i sin 0 2 sin 0 3 sin 04 . 

We choose orientation of the axis so that 

Tr(yV 7 r 7 s 7*7 u 7c) = -8e pqrstu , e Tazxyb = i. (B.2) 

With this kinematic configurations, let us write down the expression for the various 
energies which will occur in the expression for the propagators. 

E P'+P = (Ip1 J + \p\ 2 + 2 \p\\p'\ cos0i)T (B.3) 

E p i +k = (|p '| 2 + \k\ 2 + 21A;||p'| sin0i cos0 2 )^, 

E p f +v+k = (|p'( J + \k\~ + |p| 2 + 21A:||p'| sin0i cos0 2 + 2|p||p'| cos0i)T 

For convenience in notation we will also define 


iuj- = itUm + e{p — p c ), ujm = (2m + 1)7 tT. 


(B.4) 
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Now substituting the expression for the propagators from (3.21), into (4.15), we 
obtain 


d b p' 


If (P + k)f (-fc)f (-„)> = ^E/(As x (B.5) 

| ^ Tr[7 a 7 Q 7 a: 7 / 3 7 J: 7 l/ P + ]A t (ia;_,p / + k + p){p’ + k + p)(t, a )X 
\t,i;,u=d= 


\(w_y +p)(p' + i%,«)A 1 ) (A;_,.pV( i ^)J + 

(Exchange term). 

The exchange term is the second term in equation (4.15), arising due to the second 
Wick contraction. We can obtain the exchange term from the term written down by 
the performing the following operation on the first term written in (B.5). 


a —> z, 

(p +p + k) —>■ (p —p), 


z —y a, 

(p + p) —> (p — p — k). 


(B- 6 ) 


Performing the trace over the 7 matrices, we obtain 

Tr['f'y a 'f'y^'f'y v P+W + k + P)(t, a )(p' + p)(u,p)(p%,v) 
'Atuvp' a (p'+p) z (p' - q )a 


(B.7) 


Epi E(pi + p) E( p i_ q) 


+ permutations in (a, z, x ) 


4 1 u(k y p' b + p b k y ) 4 u v p' y p b 41 v(k y p' b - p b p' y ) 


E p ' +p+kEpi _|_p EpiEpi +p Epi +p+ kEpi 
The next step is to perform the following four integrals 


-M 


ipbik.y p b k y 2 3 /3 
—A 2 _ e 3 

ip b iky p b k y 2 3 (3 

-A 3 _ e 3 
ip b iky p b k y 2 3 /3 

—A 4 e 3 

ipbiky p b k y 2 3 f3 


E 


E 


E 


E 


0_ ^ ( WM + pW + P + k), + permutation ^ x 


t,v,u=± v Ep'E(p'+ p )E(pi_ ( p 

A t (iuj-,p' + k + p)A u (iuj- 1 p + p)A v (iuj-,p'), 




d 5 p ‘' 


E 


4 1 U Pb ky 


(27t) 5 . ^ ^ 


X 


t. V .u=± ^P'+P+k^P'+P 

A t {iuJ-,p' + k + p)A u (iuj_,p' + p)A v (iuj_,p'), 


E 


4 tukyP'h 4 tv k y p' b 


E r J 4-rt4-k- Er\' 


X 


d 5 p' 

(2 ttj j r u \Ep'+P+kEp'+p JZp'+p+kiZp' 

A t (iuj-,p' + k + p)A u (iu}_,p + p)A v (iu- 1 p), 


d 5 p r 

(27rp 


4 uvp' p b 4 tvpbp' 


Ep'E p i + p 


+ 


E 

t,V,U =zb 

A t (iu-,p' + k + p)]A u (iw_,p' + p)A v (i(j_,p'). 


Ep' -\-p-\- kEp' 


(B. 8 ) 
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Here we have divided by the external momenta and multiplied by a negative sign 
in each term so that we obtain terms directly related to the transport coefficient of 
interest. Each of the terms involve a sum over Fermionic Matsubara frequencies. 
They are evaluated using the following expression 

+ p + k)A u (iuj-,p' + p)A v (iuj-,p') (B.9) 

m 

_ tf^-^p'+p+k tcfjj— ) V'fi^E'p'+p we/m) 

(tE p i +p+k vE p >) (tE p i +p+k uE p ' +p ) (uE p ' +p / E p ,— p +k ) (^ E p t ^ p vE p i ) 

H _ vf(E p / - yep-) _ 

(vEpi l E p /+ p ^k ) (v E p ! uE p > +p ) 

= M. (B.10) 


In arriving at the above sum, we have neglected terms independent of the temperature 
and chemical potential. An important observation from (B.3), is that all terms in 
the Matsubara sum is independent of the angles 0 4 due to our choice of directions 
of the external momenta. Let us first evaluate A\. The numerator is proportional to 
P'aPxP'z- After substituting for these components from (B.l), we that the numerator 
is proportional to sin </> 4 cos </> 4 . Then on performing the angular integral over </> 4 , this 
term vanishes since each term in the Matsubara sum in (B.9) is independent of </> 4 . 

The term A 2 term is given by 

d 5 p\ 4 tu[pbk y 


ip b ik y p b k y 2 3 (3 


Y' / » Pl y' x ^ 

Urn J t,V,U=± E Pl+P+k E p 1+ p 


(B.ll) 


We perform the sum over a pair from {t, u, v} in each of the 3 terms in the Matsubara 
sum JH as follows. For the first term in the A4 we sum over the pair u,v, for the 
2nd term the pair t, v and for the 3rd term we sum over the pair t, u. We can then 
change variables of the p' integration so that the numerator in each of the terms 
in the Matsubara sum is a function of the integration variable //. For instance in 
the first term we substitute p' + p + k = p', while in the second term we substitute 
p' + p = p'. This ensures that the numerator of each term in A4 is a function of 
E p i with p' being the integration variable. After performing these manipulations we 
obtain 


—A 2 

ipbiky 


J2 J ^PifiEp' ~ ten-) 

1 | 1 

W V ' ~ Epi +p ) (E^, - E* +k+p ]) + (El, - El,_ k ){El, - Ep, +p ) 

i 

IE' - EEEE - WEE)' 


(B.12) 
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We can now take the p, k —> 0 limit and expand each term in the above expression as 
a series in p, k. The order of limit in these expansions does not matter. The leading 
singular term and finite terms are 


r —A 2 2e 
Inn - = -— . 

p b ,ky^o ip b iky (27r) 5 


x 


y: / d 5 pt f (E p i — tep,- 

t=± ■' 


(B. 13) 


8 p' 4 cos 2 p i sin 2 0! cos 2 p 2 4 ptp\ cos 2 pi sin pi cos < 

+ 0 (p b , k y ). 


Note that the singular term in pb on (j ) 2 —> 7r— (j ) 2 and therefore vanishes on performing 
the angular integration over </> 2 using the ie prescription. Therefore one is left with 
the following expression 


lim T —t~— 

Pbyky^O ipbiky 


2 e 3 


8(2vr) 


E 

t=± ■ 


dp d(f)id(f) 2 d(f) 2 ,d(t)A x 

sin 0 3 i/(|p'| - 


(B. 14) 


sin (pi sm 


cos 2 (pi cos 


The angular integration using the ie prescription can be performed using the following 
results 




2 . 


(B.15) 


Substituting these values for the angular integrations we obtain 


lim 

Pb,ky^> 


—A 2 
o ipbiky 


r°° e 3 

/ ^'lfor 3 IftP' ~ e V~) ~ f(p' + e V-)\ 
0 

e 4 /i_ 

167r 3 


(B.16) 


Proceeding in a similar manner, we find that when we sum over pairs in {t, u, v} as 
discussed earlier in term A% it vanishes. The term A 4 also vanishes due to the same 
reason. Now the exchange term in (B.5) also contributes equally. Therefore using 
the transport coefficient in the Euclidean theory is given by 


cT> 


3e 4 

647t 3 ^ 


(B.17) 


Note that from (4.3) the transport coefficient is 3/8 times the ( jjj) correlator. 
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B.2 Contributions to 

After performing Wick contractions on the correlator corresponding to this transport 
coefficient we obtain 

° 2 r d 5 p' 


<J°(P + k)f(-k)T”(-p)) = i— Y, 


2 *P^-J ( 2 nf 


E Tr ([7“7“7“'7' 5 (7 T ip« + iu-Y)l V P+] X 


t,v,U= zb 


+ p + k)(p' + p + k) {oiit) A u (iu_,p' +p)(p' + p) { p, u )A v {iu_,p')p' , 

d 5 p’ 


E (M7V« + ^-7*)7“tV7V-P + ]x 

m ' ' t,v,u =dz 

Attiw-.p'ApJfp' -p)(„,t)A„(iw-,p' -p- fc)(p' -P-*:)(,),„)A„(iw_,pV(«,»)) • 


(B-18) 


The second term in the curved brackets is obtained by the the exchange Wick con¬ 
traction. Examining the terms we see that we need to evaluate the following integrals 

” 2 r d 5 p ' 


E / (|y E 'tt[7“7“7V7V-P + K X 


(B.19) 


A t (iu-,p' +p + k)(p' +p + fc) (Qi i)A u (io;_,p' +p)(p 7 +p) (j 8,u)A 1 ,(ia;_,p , )p , ( V ,«)> 

d 5 pi 


Ao = i 


^E/p E ± Tr W^Y^]®- x 


A t (iuj-,p +p + k)(p +p + k)( ait )A u (iu-,p' +p){p + p)u3,u)A v (iuJ-,p')p{ UtV y 

Evaluating the trace over the 7 matrices, we end up with 


• e /* d p' 

I_ ’wE/ (W,.? 


d 5 p' ^ 4 tuvkyPb p? 

( 2 71 ") 5 E p iEpi+p+kEp'+p 


x 


(B.20) 


A t (iu}_,p' +p + k)A u (iu_,p' + p)A v (iuj_,p'), 

d 5 p' ^ (ttu(k y p' b + p b k y ) Auvp'yPb 


. e 2 y^ f d p' 

2 ( 2^55 E 


(2 tt) 5 twu=± \ Ep/_|_p_|_fcEp/_)_p EpiEpi +p 


Atv(kyp u - p b p'y) 


iu)-A t (iu-,p' +p + k)A u (iu-,p + p)A v (iu-,p). 


E p ' +p+kEp' 

The sum over Matsubara frequencies in Ai is performed using (B.9). To perform the 
sum in A 2 , we use the following 

1 y^Juj-A t (iuj-,p' + p + k)A u (ioj-,p + p)A v (ioj_,p) 


P 


Ep' +p +/,‘f ( E p i^p +k tep ) 


+ 


Ep'+pf ( E p i +p — uep,-) 


(tEp'+p+h v Ep'^) (tEp'+p -)_/- uEp/j^p) (uEpt+p tEpf-^p+k^uEp'+p vEp/^j 
Ep/f (E p f - veil-) 


{vEpt t Epi-^.p-^f^j {v Ep/ uEp/^p^j 


(B. 21 ) 
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We can now proceed as discussed in B.l to perform the sum over {t,u,v} and then 
change variable of integrations. Finally after performing the angular integrals using 
the ie prescription, we obtain 


—A\ . e 

inn 7 — : — = % 

Pb,ky—^0 ip h lky 

-A 2 


32tt 3 
e 


lim -— = i- 

Pb,k y ^o ipbiky 327T 3 


pOO 

Y / dp' p'f (p' — ten~) 

t=± Jo 

poo 

Y / dp pp f(p - ten-). 
t=± J o 


(B.22) 


The exchange term in (B.18) also contributes an equal amount to the correlator. 
Therefore, we have 


lim 

Pb,ky^ 0 


j a (k + p)j x (-k)T tz (-p) Ai + A 2 

-- = inn —2 


ipbik y 


Pb,k y ^fO ipbiky 

e 2 r°° 

= i 8^3 Y J d PPf(P ~ te P~)i 

7 r 2 T 2 


= 1 


167T 3 


(e 2 (/i -) 2 + —)• 


(B.23) 


.( 6 ) 


B.3 Contributions to Q 

The correlator corresponding to the transport coefficient also involves contact 
terms. We will first evaluate the contribution from Wick contractions of the flat space 
currents in the correlator . We then proceed to evaluate the contribution 

from the contact term. 


(j a fl ( P +k)Tjn-k)Tjn-p)) 

e ^ /" d 5 p' 


(B.24) 


32/3 


Y / 7^)5 Y ^b a i a {n T w' x + i x ^-)i p {n T W'z + i z ^-)YP+] 

™ ' ' t,v,u=± 


x(p' +P + k){t, a ){p' +p)( u ^)p{ VjU) A t (iuj_,p' + k + p)A u (iuj_,p' + p)A v (iuj_,p') 
+ (Exchange term). 
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Among the four possible traces it can be seen that only 3 of them contribute, they 
are given by 


Ai — 


re 


E / 7535 E Tr[ 7 V7V7V<P+]p>- 

m n m -I- 


32/3 


A 9 — 


E / E 

m ' ' f it 1/-I- 


A* — 


t,v,u ==b 

x(p' +P + fc) ( i,a)(p' + fc + p)A u (wn_,p' + p)A„(mu_,p') 

+ (Exchange term), 
ie ^ f d 5 p ’ 

32/3 „ 

m v y 

x(p' + p + k)(t,a)(fi +P)(u,/3)P\v,v)A t (iuJ-,p' + k +p)A u (iuj_,p' + p)A v (iu_,p') 

+ (Exchange term), 

A)E [/% E Tv^^yv-p+i^-) 2 

m ^ y fc,v,n=di 

x(p' +P + k)(t,a)ip + P)(u,/3)P{ v ,v)A t (iu-,p' + k +p)A u (iuj_,p’ + p)A v (iuj_,p’) 

+ (Exchange term). 

(B.25) 


Now we follow the methods adopted to evaluate the previous correlators. We first 
perform the traces in each of them and follow that up with the Matsubara sum and 
then take the limit, pb, k y —y 0. Lastly we perform the angular integrals with ie trick. 
We find that the terms Ai and A 2 contribute equally and they are given by 


t4i 

a 2 


E f d 5 p 1 4 tuv k v p b Pr iu- ... , . . . . , . , .. 

jAi E F F r - A,(p' + fc+p)A„(p'+p)A„(p') 

J ( tvu= ± Ap/ L/pi+p+k-E'p'+p 


32/3 

+ (Exchange term), 

E 


d 5 p' 4 tuv k v p b p'z iuj_ . . , , . . . . . . . .. 

Y - - yFb % - A t (p'+ k+p)A u (p'+ p)A v (p') 


—e 

t,v~u^=± Ep'Ep’+P+kEp’+p 

+ (Exchanged term). 


(B.26) 


The Matsubara sums are performed using equation (B.21). After a change of vari¬ 
ables of integration and taking the external momenta to zero, we obtain 


lim 7 ——^- 

Pb,ky—>0 ipfoifcy 


2 e 

647T 3 


poo 

Y / dp'tp , 2 f{p' -ten A). 
t=± J 0 


(B.27) 


Here again we have used the ie prescription to evaluate the integrals. The factor 
of 2 results from taking into account the equal contribution of the exchange term. 
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Evaluating the A 3 term involves performing the following Matsubara sum 


1 

1 


Y( iuj -) 2 A t (iuj-,p +p + k)A u (iuj-,p + p)A v (iu-,pi) 

m 


_ ^^p'+p+kf (-^p'+p+k ten A) u Ep'+pf {Ep'+p uep,-) 

(tE p / +p+k vEpi ) (tE p /^_p_^ k uEp' +p ) (uE p t _(_p f-^p'+p+fc) i^Ep' +p vE p f ) 

_ vE 2 ,f(E p i - yen-) _ 

(vE p i tEp^p^A'j (vEpt uE p > +p ) 

= AT. (B.28) 


Performing the trace over the 7 matrices and substituting the Matsubara sum A/", 
we obtain 


e y- f d 5 p ’ (4:tu(kyp , b +p b ky) _4:Uvp' y p b 

32 /3 m J ( 2?r ) 5 t,v^± ' E p f +p+k Epi +p EpiE p i +p 

Atv (ky p’b — PbP'y) \ 

- F - F — X U 

+ (Exchange term). 


(B.29) 


Now summing over {t, u, u}, making appropriate shifts in the integrations and taking 
the zero momentum limit, we get 


~A:> 


lim 

p b ,k y ^o ip b ik y 


2 e 

647T 3 


poo 

Y / dp tp 2 f{p 
t=± J o 


te p-). 


(B.30) 


Therefore the putting terms A 1 , A 2 , A 3 together we see that the contributions from 
the Wick contractions of the flat space currents is 


lim 

PbPy^O 


{3 a fi(-p- k ) T fi{k) T fi{p)) _ 3e 


ipAky 


32tt 3 


poo 

Y / dp p r2 tftp’ -ten-). 

t=± 


(B.31) 


We will now obtain the contribution from the contact terms. For this we first 
Fourier transform the interaction S 2 E given in (4.11), this yields 


5S 


( 2 ) 


dh XT (k)Sh ZtT (p) 


16/3 ^. 


d 5 p' 
(27r) 5 


fk,p+ k + p)tPbi1 X T + k y ~i y Yl x )P-f){u m ,p'). 

(B.32) 


Considering the Wick contraction of the above contact term along with the flat space 
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charge current, results in the following two terms 


—~B = Bi + B 2 

XPlyXhy 

£ti = —i y 

16 k y (5 

* m 


(B.33) 


B 2 = 


d b p' 

(2 

d 5 p r 


Tr {rS(Lu m ,p + k + 3 /) 7 i 7 *f%,y)P + ], 


IQpbP 


TY[rS(bu m ,p + k + p') 1 b 1 z l x S{u m ,p')P + } 


We proceed to evaluate i?i. Substituting the expression for thermal propagators we 
obtain, 

lim lim B\ — lim V TrhVbVbV^+l 

Pb,ky->OpO,ko-tO Pb,ky-^0 ldky /3 J (2l l ) 5 

^ m. x 7 /.?/—+ 


X T (p' + p + k) (t,a)P' (u,/3) A * («*>- , p' + p + k) A u (ioj_ ,p'), 


= lim — 

Pb*ky —^0 16 ^ 

t,U=± ' 


jV t tf(E p , +p+k - te/i_) - uf (E p i - ue/i_) 

( 27 t ) 5 uEpi 


= j dp> rf< M < M<ft 3 (27r) ^ sin 3 0 i sin 2 0 2 sin 0 3 

x - 7 “—— + . „ 1 2 , ) + 0(k y /p b ,p b ). (B.34) 

\ cos 4p /2 cos 2 (pi J 

To obtain the last equality we have summed over {t, u } and then taken the external 
momenta to zero. We have taken k y — > 0 first and then taken the limit p b —» 0. But 
the final result for the non-singular term is independent of the order of limits. Also 
note that the singular term in p b is odd under (pi —M t — (pi. Therefore it vanishes on 
performing the angular integration over (pi. If the order of limit is interchanged the 
resulting singular terms in k y also vanish on angular integration and the final result 
for the non-singular term is invariant. The angular integration in this term can be 
evaluated easily by the ie prescription, 


r sin 3 0 i 
/ 0 COS 2 01 


= -4. 


(B.35) 


Substituting this result into £>i we obtain 

p r°° 

lim Si = dp tp' 2 f{p' — tefi-). (B.36) 

Pb,ky^ 0 647T J J o 

Going through a similar analysis on B 2 it can be shown that the result for B 2 is 
equal to that of B { . Thus the total contribution from the contact term is given by 

dp tp' 2 f {p — tep-). (B.37) 
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Therefore putting all the terms from A and B together we obtain the final result for 

the correlator of interest to be 

—1 —e . r°° 

^ti“(p+fc)r™(-*or"(-p)) = J d P 'tp a f( Pl -ten _) 

= ^(eV-WUV). (B-38) 

B.4 Contributions to A 3 6 * 

There are two contributions to this correlator, one from the direct Wick contractions 
of the flat space currents and the other from the contact terms. Let us first evaluate 
the former contribution, performing Wick contractions leads to 

(TJt(p + k)T;n-k)Tjt(-p)) = jElJ^ x ( R39 > 

(Tr [7 T ip' a + j a iu_)^r a (Yip' x + 'fiu>-)'yP('y T ip' z + 

+ p + k)A u (iuj_,p' + p)A v (iu_,p')(p' + p + %, a) (p' + p){u,p)P[ v , u )) 

+ (exchange term). 

Expanding the above express it can be seen that there are 4 terms which are non¬ 
vanishing these are given by 

x 

m ' ' 

A t (iu-,p' + p + k)A u (iuj_,p + p)A v (iw-,p')(p' + p + k)(t, a )ip' + P)(u,0)p{ v „) 
+ (exchange term), 

/ |^Tr[ T vv/T Ywy x 

A +p + k)A u (icj_,p' + p)A v (iuj_,p')(p' + p + k)( t>a )(f + p)( u ,0)P'( v ,u) 

+ (exchange term), 

^ YW- X 

m ' ' 

A t (iu-,p' + p + k)A u (iuj_,p + p)A v (iu-,p){p + p + k)( t , a )(p' + P)(u,p)P'( v , v ) 
+(exchange term), 

A ‘ = -659E / I^tVtVtVJ’h-K- 

A t {icj_,p' +p + k)A u (iu_,p' +p)A v (iuj_,p')(p' + p + k) t ^(p' + p) u ,pP f ( v ,^ 

+ (exchange term). (B.40) 

We find that contribution of the terms A 1 ,A 2 and A 3 are equal. Furthermore the 
contribution of the direct Wick contraction is equal to that of the exchange term in 
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each Ai, A 2 , A 3 , A 4 . After performing the trace over 7 matrices, we find that Ai is 
given by 


Ai 


2 i y-~^ f d 5 p' y-~^ Atuvp'a 2 

64 P m J (2vr) 5 t u ^ =± Ep'+p+pEp' +p E p r 


A t (iu-,p' +p + k)A u (iu-,p' + p)A v (iu>-,p'). 


(B.41) 


ffere the factor of 2 is due to the contribution 0 
in discussed in the previous examples we obtain 


lim — -—— 

Pb,ky^>-0 ipbiky 


2 i 
’128tt 3 


r°° 

Y / dp'P 3 f(P — te H-). 

-- - t=± J o 

Since the contribution of A 2 and A% is identical to that of A\ we have 

A 1 A 1 A a A _ poo 

Y / d P'P >3 f(p' -ten-). 

t=± Jo 


Al + A 2 + A3 KJU 

hm-=-- 

p b ,k y ^o ipbiky 1287T 3 


(B.42) 


(B.43) 


To evaluate the term A 4 we need to perform the matsubara sum which is given by 
4 Y( iu -) 3 A t( iu -’p' + P + k)A u (iuj_,p' + p)A v (iaj_,p') 

“ m 

Ep'+p+kf (Ep'+p+k — tefi-) ^ Ep, +p f (Ep/+p — uefi-) 

vEpi') (t-E'p'+p+Zc ti-Ep'-j-p) {uEpi+p (uEpf+p vEp> ) 


(tE’pZ+p+fc U J-Jp‘ 

_ E p'H e p' ~ ve V~) _ 

{vEjpf ^-Ep'+p+fc) (vEpf , n£'p/_|_p) 

At the end of the usual manipulations, we obtain 

lim = -AA [ dp'p' 3 f(p' -te^). 

Pb,ky ^0 1287T 5 ^ Jo 

Therefore the contribution from the Wick contractions of the flat space current is 
given by 


(B.44) 


(B.45) 


hm 

Pb,ky^ 0 


(■ T Ta (p + k)T TX (-k)T TZ (-p)) 
ipbik y 


-8 i 
128vr 3 


Y f dpp i f{p 

t=± d 


t e /i_). 

(B.46) 


Now lets examine the contact terms. As seen in the equation (4.37), there are 3 
contact terms. We will evaluate one of them, it can be seen that the contributions 
from all the 3 contact terms are equal. Performing the Wick contraction of the 
contact term with the flat space stress tensor we obtain 


—Bi 

ipbiky 


“3^54 / (^ Tr ((7 Ti Pl + + p + f)(7Vf )5(w m ,p')P + ] 

+—Tr[(7 T ^l + 'fuo-)S{u m ,if + p + k)('y y Yl x )S(u m ,p , )P + \) ■ 

Pb J 

(B.47) 
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We then follow the same steps as before to evaluate this integral and get 


B i 


lim 

Pb^-y—^O iPblKy 


-8 i 


3 128vr 3 


poo 

Y / dp j/ s f {'p — ten ). 

t=± 


(B.48) 


The other two contact terms B 2 and B 3 also give the same amount of contribution. 
Therefore the total contribution from the contact terms is given by 


B 


8 i 

128vr 3 



dp p 3 f{p' 


t e /i_). 


(B.49) 


Combing the terms A and B we see that the value for the three point function of 
interest is given by 


, ( T Ta (k + p)T TX (-k)T TZ (-p )) 

Inn-;— - 

Pb,ky—>0 IPblKy 


_i aj _ r°° 

fluffs Y J 0 dpp 3 f{p-tey _). (B.50) 
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